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Interfaces pinned by quenched disorder are often used to model jerky self-organized critical motion. We 
study static avalanches, or shocks, defined here as jumps between distinct global minima upon changing an 
external field. We show how the full statistics of these jumps is encoded in the functional-renormalization-group 
fixed-point functions. This allows us to obtain the size distribution P{S) of static avalanches in an expansion 
in the internal dimension d of the interface. Near and above d = A this yields the mean-field distribution 
P{S) ~ 5'-3/2g-s/4S„ ^jigjg 5^ is a large-scale cutoff, in some cases calculable. Resumming all 1-loop 

contributions, we find P{S) ~ S~'^ eyip{c{S / SmY'^ - ^{S/Smf^ where B, C, 5, r are obtained to first 

order in e = 4 — d. Our result is consistent to 0(e) with the relation t = tq := 2 — where C, is the static 
roughness exponent, often conjectured to hold at depinning. Our calculation applies to all static universality 
classes, including random-bond, random-field and random-periodic disorder. Extended to long-range elastic 
systems, it yields a different size distribution for the case of contact-line elasticity, with an exponent compatible 
with T = 2 — to 0{e = 2 — d). We discuss consequences for avalanches at depinning and for sandpile 
models, relations to Burgers turbulence and the possibility that the relation t = r^he. violated to higher loop 
order. Finally, we show that the avalanche-size distribution on a hyper-plane of co-dimension one is in mean- 
field (valid close to and above d = 4) given by P{S) ~ Ki {S)/S, where K is the Bessel-ZC function, thus 

'Hiypcr plane — g • 



I. INTRODUCTION 

A hallmark of complex non-linear systems, as well as sys- 
tems with quenched disorder or inhomogeneities, is that the 
response to an applied field is very often not smooth but in- 
volves jumps, bursts or avalanches. This is true for domain 
walls in a magnet responding to a change in external mag- 
netic field, leading to the Barkhausen noise lU |2l, the flux 
lattice in type-II superconductors upon varying the field |[3| 
the contact line of a liquid partially wetting a disordered sub- 
strate when emptying the container 10] [5], or piles of granu- 
lar material when adding grains Q. Other examples are mo- 
tion of cracks in brittle materials, dry friction |7| and earth- 
quakes fSl 19). This jerky behaviour often arises as a non- 
equilibrium phenomenon, but it may also occur at equilib- 
rium, in systems with many metastable states, as they switch 
from one global minimum to another when an external pertur- 
bation is applied. The statistics of these jumps ubiquitously 
exhibits scale invariance and self-organized criticality as in 
sandpilesi 10|, with power-law tails for the probability of rare 
large events: if one defines an event size S, the probability 
distribution behaves as P{S) ^ S~'^, up to some large-size 
cutoff Sm, e.g. imposed by the finite system size. 

An outstanding question is the degree of universality of the 
jump statistics. This can in principle be answered in partic- 
ular prototype models, but even then it turns out to be rather 
difficult to obtain analytical results. One such class of mod- 
els are sandpile automata which are dynamical systems where 
small events can trigger large avalanches ifTOl [TTl [121 [T3l . 
There, some beautiful results have been obtained analytically, 
but the full avalanche statistics, including the distribution of 
avalanche sizes P{S), has not yet been obtained in the cases of 
physical interest, such as spatial dimensions d = 2,3. Mean- 
field theories |14, 15| predict r = 3/2 and various scaling 
arguments have been constructed lT0l[TTl[T2l[T3l[T6l[r7l . not 
all in mutual, agreement especially for d = 2. Among these. 



the conjecture t = 2 — 2/d seems to be the best guess lfT8l[T9]| 
in d = 3, but most often one has resorted to extensive numer- 
ics. 

Another class of prototype models are random field Ising 
magnets (RFIM) in an external field H. There one may study 
either the changes in the ground state as H is varied or the 
non-equilibrium zero-temperature dynamical evolution from 
an initial state. Tuning the system near the critical values of 
field and disorder where the macroscopic magnetization jump 
vanishes, one can study finite-size avalanches, and these are 
found to exhibit self-organized criticality. Extensive numeri- 
cal work has been caiiied out to determine avalanche-cluster 
statistics |20|. In parallel, the avalanche size distribution was 
studied [21 1 using a a field theoretic RG approach in an ex- 
pansion in d = 6 — e. However, since the RG used there is 
based on the dimensional reduction property, which is well 
known to fail for even simpler disordered problems, as well 
as for the static RFIM, the status of the result of Ref. 121], 
T = 3/2 + O(e^) remains to be clarified. Furthermore it 
has been recently argued, mostly on the basis of numerics, 
that the avalanche statistics for the RFIM ground state and for 
its non-equilibrium dynamics are not distinct but belong to 
the same universality class [22 1. Hence, another outstanding 
question is to clarify the possible differences between equilib- 
rium and non-equilibrium avalanches in the RFIM, and in a 
broader class of models. 

Elastic systems pinned by a quenched random substrate 
provide yet another frequently used model for avalanche phe- 
nomena. In the statics they are known to exhibit glassy phases 
where the Gibbs measure, while localized near a minimal en- 
ergy configuration, does exhibit jumps at low temperature due 
to the presence of many low-lying metastable states. As an 
external force / is applied, the system exhibits a depinning 
transition: at T = it starts moving at non-zero velocity only 
for f > fc- The motion near the depinning threshold fc has 
been studied extensively numerically and is known to proceed 
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by avalanches 11231 1241 125| 1261 . These are actually not ob- 
vious to define non-ambiguously in the moving phase, hence 
some of the literature on the subject is qualitative. Under- 
standing these avalanche statistics is however a challenging 
question, since elastic plates driven by springs exhibit stick- 
slip motion, which provide a first step to model complex sys- 
tems such as earthquakes In Ref. |8| a very simple model 
for avalanches, mean-field in inspiration, was proposed and 
solved, and yields again the value r — 3/2. Surprisingly, 
there is to this day no first-principle derivation of even this 
simple mean-field behaviour, within this class of models. Be- 
sides extensive numerics, the main theoretical result is a con- 
jecture for the size exponent near depinning, equivalent to 
T = = 2 — 2/((i+C), based on a scaling argument and some 
unproven assumptions Il27ll28l . Thus there is the need for an 
analytic tool to approach this problem. The Functional RG 
(FRG), a powerful field-theoretic method to deal with disor- 
dered elastic systems in the statics 1291 l30l [311 l32l [33l l34l [35l 
[36, 37, 38 , 39, 40, 41 , 42, 43, 44| and in the driven dynamics 
E21 [45] [S [13 (48, 49 50, 51 52, 53], as well as with random 
field models fj^i ,55j |56i i57j |58j |59l [60] [SH, has been quite 
succesful at computing for instance the roughness properties 
of the displacement field u{x), in a dimensional d = 4 — e 
expansion. Following the procedure proposed in Il62ll63l[64l . 
the numerical determination of the predicted FRG fixed-point 
functions in dimensions d = 1,2,3 have confirmed the theory 
to high precision, up to two-loop accuracy both in the statics 
|[65l and at quasi-static depinning ll66ll . This was an important 
test, given that the field theory is quite non-standard: the cou- 
pling constant is a function of the field, A{u) = —R"{u), the 
renormalized second cumulant of the pinning force, and the 
zero-temperature effective action is non-analytic, e.g. A{u) 
exhibits a cusp at m = 0. However, until now no calculation of 
avalanche distributions has been even attempted using FRG. 
In fact, possible difficulties in handling the fast-jump motion 
within the field theory have been emphasized [46 j , and it was 
unclear whether it was possible at all. 

The aim of this paper is to show that one can extract the 
avalanche statistics from the FRG in a controlled way. It is a 
priori quite involved as it requires the calculation of the non- 
analytic part (the cusp) of all cumulants of the renormalized 
disorder, which encode the jump distribution. However, it 
turns out to be feasible in the end, as we find that remark- 
ably simple self-consistent equations are obeyed by suitable 
generating functions. As an application we compute the dis- 
tribution of avalanche sizes from first principles in an expan- 
sion in d = 4 — e. We derive the mean-field result and obtain 
a 1-loop, i.e. 0(e) prediction for P{S). Here we study the 
case of static avalanches, also called shocks, which are de- 
fined as jumps between distinct global minima of the energy 
upon changing an external field. The most convenient setting 
is to add an external harmonic well with a variable central 
position, i.e. consider an elastic interface in a random poten- 
tial tied to a harmonic spring. As a function of the spring 
position, the center of mass of the interface changes in dis- 
crete jumps. We study all static universality classes including 
random-bond, random-field and random-periodic disorder, as 
well as long-range elastic systems. We obtain not only the 



exponent r but the full scaling function of the avalanche dis- 
tribution. The latter is (almost fully) universal with respect 
to small scales, and the dependence in the large scales, i.e. in 
S/Sm, is calculable. A short summary of some of our main 
results, together with a parallel numerical study, has recently 
appeared ll26l and here we provide the necessary details of the 
approach and of the calculation. Although we do not directly 
study depinning here, there are connections between static and 
dynamical FRG, the main idea of the method being similar 
As found in a companion study |67| of avalanches at depin- 
ning, the results of the present paper are most likely to hold 
also for depinning to one loop. Since for some observables, 
such as the roughness exponent C,, differences between statics 
and depinning within the FRG appear only at two loop |32|, 
the question of the difference between non-equilibrium and 
static avalanches remains open and requires a 2-loop study 
||67| . Note that to our knowledge, while avalanches near de- 
pinning have received a lot of attention, the static avalanches 
and shocks defined here have not. Their connection to the 
Burgers equation 1621 [63l . which provides their d = limit, 
helped formulate the problem in a way amenable to FRG cal- 
culations. 

Let us also point out that there are interesting relations be- 
tween the depinning in the random periodic class, i.e. charge 
density waves (CDW), and sandpile automata l24l [68l . Sand- 
pile automata are also related to spanning trees and loop 
erased random walks |[l8l[l9l[69l[70l[2ll[2Il[23[74J, and this 
has lead to recent FRG predictions for the latter l75l . Hence 
the present work is also relevant to sandpile models. An im- 
portant issue is the status of the conjecture r = t^^, which in 
sandpile models reads r — r<^=o consistent with a value ( = 
for CDW depinning. Of course, in the present paper the ques- 
tion becomes whether such a relation also exists in the statics. 
Finally, let us note that our results extend beyond disordered 
systems. Indeed in some complex dynamical systems, the ran- 
domness lies in the initial condition, and avalanches or shocks 
appear in the course of a non-linear deterministic evolution. 
This is the case for the sandpile automata (e.g. in the con- 
served energy ensemble) as well as for the problem of decay- 
ing Burgers turbulence. The latter is indeed intimately con- 
nected to the pinning problem, and to the FRG, as is recalled 
below. Therefore the methods introduced in this article should 
find applications in a much broader class of non-linear com- 
plex systems, possibly including turbulence and spin glasses. 

The outline of the article is as follows. In Section [HI we 
define the model of an elastic manifold in a random poten- 



tial and a harmonic well. In Section III we give the expected 
scaling form for shock ditributions. In Sectio n |IV| we detail 
the connection to FRG functions. In Section Mwe present 
the tree-level calculation which leads to the mean-field result. 



In Sections VI 


and 


VII 


and in Section 


VIII 


we ( 



VII we present the 1-loop resummation 



tensions to non-local elasticity are given in Section |X] More 
technical material is relegated to various appendices. 

Note that we will not review manifolds in random poten- 
tials, for which we refer to l76ir77l . nor static FRG methods 
for which we refer to fTSl for a pedagogical introduction and 
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to BSia for details of the FRG calculations. 



and for a uniform Wx = w one defines the ground state energy 
of the system per unit volume: 



II. MODEL 

We consider an elastic interface parameterized by a one 
component (N = 1) displacement field noted u{x) = Ux- It 
is subjected to a random potential V{u, x) and to a harmonic 
well centered at Ux = Wx- It is defined by the Hamiltonian: 



9q^\Uq- Wq]"^ + / V{Ux,x) (1) 



where we denote — J d'^x and the Fourier transform Uq — 

qI-q^Ux and = / -(^^fya-- A cutoff at small scale a; ~ a or 
q ~ A, is implicit everywhere. In most of the paper we focus 
on the choice of a standard local elasticity: 



9n 



(2) 



and often set if = 1 for convenience. The bare disorder is 
assumed to be short range in internal space, statistically trans- 
lationnaly invariant, with a bare second cumulant: 



V{u,x)V{u',0) ^ 5{x- x')Rq{u-u') 



(3) 



whose precise form is unimportant, apart from some global 
features which determine the universality classes, mainly (i) 
random bond (RB) with Ro{u) a short range (SR) function, 
(ii) random field (RF) with Rq{u) ~ — a long range (LR) 
function, (iii) random periodic (RP), with Ro{u) a periodic 
function of period (arbitrarily) set to unity. 

To study the equilibrium statics we define the renormalized 
potential V as the free energy of the system, and the renor- 
malized force as its functional derivative w.rt Wr,'- 



-Vlw]/T _ 



DUr& 



-■H[u-w]/T 



yx[w] = I dxx'i^x' - K')) 



(4) 
(5) 
(6) 



where (. . .) denotes thermal averages over 7i[u; w] in a given 
disorder realization and [. . .] is reserved for arguments of 
functionals, while (...) is used for arguments of functions. 
For further details and notations we refer to ||62| 1631 1641 |65l 
[66 i where this model and the observables where introduced 
and studied. 

In this article we are mostly interested in energy minimiza- 
tion as w is varied in a given realization of the random poten- 
tial V{x, u), i.e. the T ~ problem. At T = the minimum- 
energy configuration is denoted and Ux{w) for w uni- 
form, thus dropping the expectation value. V[w] becomes the 
minimum energy. 



V{w) := L-'^Viiwx = w}] 
L'^ being the volume of the system. Its derivative w.r.t w: 

V' {w) = m?'{w — u{w)) 
u{w) = L-'^ / Ux{w) 



V[w] = min?i[u; w] 



(7) 



(8) 

(9) 
(10) 



coincides with the force per unit volume exerted by the spring. 
We are ultimately interested in the small-mass limit m ^ 
where scale invariance becomes manifest. As shown previ- 
ously ||30l|65ll26l|72l|79l|80l, the optimal interface is statisti- 
cally self-affine with {u{x) — u(0))^ ^ jxp^ and a roughness 
exponent ( which depends on the class of disorder, and with a 
e = A- d expansion |32|: C = e/3 for RF, C = for RP, and 
C = 0.2083e + 0.00686e2 for RB (C = 2/3 in d = 1). This 
holds for scales Lc < L < where Lc is the Larkin length 
(here of the order of the microscopic cutoff) and ^ 1/m, 
the large scale cutoff induced by the harmonic well. It is use- 
ful to picture the interface as a collection of {L/LmY^ regions 
pinned almost independently. 

Note that in this article we study the static problem where 
the interface finds the global energy minimum for each w. The 
function u{w) is then a single valued monotonically increas- 
ing function of w. One can also study the quasi-static dynamic 
problem where w{t) grows very slowly as a function of time, 
and the interface visits a deterministic sequence of metastable 
states 1641. In that case u{w) is history dependent, although 
due to the no-passing theorem |81 1 different initial conditions 
converge to the same asymptotic trajectory as w is increased 
by an amount larger than |66|. In the to limit this 
is a method to study avalanches near the depinning transition. 
It is different from the more standard method where the force 
is increased infinitesimally while remaining below threshold 
/ < /c, and in a sense it is a cleaner method since it pro- 
duces a steady state for avalanches without getting closer to 
threshold, i.e. without changing the cutoff length which re- 
mains i„i ^ 1/to. The study of this case requires dynamical 
FRG and will be performed in 1.67 J . 



III. OBSERVABLES 

We now define a few useful observables and notations, fo- 
cusing on the simplest case of a parabola centered around a 
uniform Wx = w. We also discuss their expected scaling form 
in the limit m 0, confirmed later via the FRG analysis. 



A. Energy fluctuations 

The renormalized potential defined in Q is a random func- 
tion of w and of the bare random potential V{u,x). Hence 
one can define its cumulants as averages over the bare random 
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potential: 



V{wi)V{w2) = L-^Riwi - W2) 



(11) 



We denote S'*^"''(it;i, . . . , w„) the same expectation values 
for non-connected averages and use that R — R. Note 
that connected averages scale with system size as indicated, 
while non-connected do not. In the limit m ^ one ex- 
pects that they take the scaling form S'^"''(wi, . . . , w„) ~ 
^d-«e^(«)(^Cy^^^ . . . , m^Wn) with 6 = d-2 + 2C. 

B. Force fluctuations 

One defines the cumulants of the renormalized pinning 
force: 



m {u{wi) - wi){u{w2) ~ W2) 



L~'^A{wi - W2) 



m - Wi) . . . {u{Wn) - Wn) 

- {-irv'{w,),...,v'{wS 

^L-("-i)'^(-l)"C(")K,...,u;„) 



(12) 



(13) 



and similarly with C*^"^ for non-connected averages. (The 
same remark as above applies for the system size depen- 
dence). In the first line of Eq. (12 1 the connected av- 
erage can be replaced by a non-connected average, since 
u{wi) — wi — (as shown by parity and statistical trans- 
lational invariance of the disorder). One has the relation 
C'(")K,...,«;„) = ...a,,„^(")K,...,u;„) 
and A{w) — —R"{w). In the limit m — > one ex- 
pects that they take the scaling form C''"-'(z/;i, . . . ,Wn) ^ 

^2n-(Ti-l)(i-nC(;(n)(^Cy;^^ ^ m''Wn), and 



(14) 



where a convenient choice of the constant Ad is given below, 
see (|B3ll. 



C. Shock observables 

Consider a uniform w. It is reasonable to assume (and 
confirmed by numerical studies f26', '65]) that Ux{w) consists 
of smooth parts, which become constant in the scaling limit 
Ux{w) ^ mr'^ as TO — > 0, and jumps, also called shocks, or 
static avalanches, that it can be decomposed as 



Ux{w) 



sfe{w - w,) 



S,^ I Sf 



(15) 
(16) 
(17) 



Here 5*^ is the "size" of the shock labelled i, and 0{x) the 
unit-step function. To each environment corresponds a unique 
set of Uj.{w), and a unique set {{wi^Si)}. From this one 
defines the normalized probability of an (infinite) sequence 
. Wn, Sn', . . .). One can also define the 1- 
point probability-density 



p{S, w) = ^ S{S - Si)S{w - Wi) , 



(18) 



so that the total number of shocks in an interval of size w is 
iV^ = / els' p{S, w). We will assume that A^^ is propor- 
tional to w. This is equivalent to 



p{S,w) = pQPiS) = p{S) 



f 



dSp{S) (19) 



where P{S)<:\S is the normalized probability that a given 
shock has a size in the interval [S, S+dS] and p^dw is the av- 
erage total number of shocks in an interval dw, i.e. the shock 
density, assumed here to be a finite number Note that the mass 
TO provides a scale ensuring convergence at large S, while the 
notation S'min refers to the avalanche size at the small scale 
cutoff of the model, which in some cases will be explicitly 
needed to ensure convergence at small S of our scaling forms, 
see below. The statistical translational invariance of the disor- 
der, together with parity implies: 

u'{w) EE -jj^ S.,S{w -Wi) = I dSSp{S) = 1 . 

(20) 



This also gives 



L-'po{S) = 1 , 



(21) 



where here and below (5") := /d55"P(5) are normalized 
moments. 

We now consider the m ^ limit. One expects, and later 
verifies that the shock distribution takes the following scaling 
form (see Il24ll26l^: 



p{S) = L'^mPS-^p{Sm'^+<) 



(22) 



where p{s = 0) is a constant and p{s) has a fast decay to zero 
at large s, suppressing large avalanches. This form involves a 
priori two exponents p and r. However if one assumes r < 2, 



which is found to always hold here, the constraint (20 1 implies 



a relation between the two exponents obtained by writing, as 

TO — !■ 0: 

/ dSSp{S) = to''-(2-^)('^+C) / dsss-^p{s) = 1 



p=(2-r)(d+C) 



/ dss^-^p{s) = 1 , (23) 
Jo 



It also implies a constraint on the scaling function. From now 
on we denote Sm the scale at which the avalanche sizes are 
cut by the mass m. It behaves at small mass as 



= cm-(^+0 



(24) 
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For now it is defined up to a multiplicative constant c, for 
which we make a convenient choice below. 

We must now distinguish two cases where the distribution 
is qualitatively different: 

(i) r < 1 : 

In that case the total shock density 



po = L m 



(d+O 



dss ^p(s) 



(25) 



is given by a convergent integral at small s — cS/ Sm (it al- 
ways converges at large s). Hence the normalized avalanche 
size distribution P{S) takes the simple scaling form 1891 : 



P{S) = S-,'p{S/S^ 
p{x) 



'^x p{cx) 



(26) 
(27) 



controlled by a unique scale S„ 
for the probability satisfies: 



The scaling function p{x) 



p{x)dx — 1 



xp{x)di 



1 



(28) 

It is itself a probability (therefore the notation p{x)). 

Two known examples with r < 1 are (for review see 1631 ): 
(i) the Sinai model, which corresponds to the d — limit of 
the random field RF case, i.e. ( — 4/3, and also to the de- 
caying Burgers equation with initial uncorrected velocities. 
There the scaling function p{s) is known, with r = 1/2 and 
p{s) ^ s^/"^ exp(— as-^) at large s (ii) d = short range ran- 
dom potential, also known as Kida turbulence in Burgers lit- 
erature, which corresponds to ^ = 1. There p(s) = \se' 
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(ii) 2 > r > 1: 

Here the total shock density in the small-mass limit S-^ 
Sm is controlled by small shocks: 



Po 



p{S)dS 



(29) 



The normalized size-distribution can then be written for S ^ 

•^min is: 



P{S) = C{S)S-,\S/S,nr^f{S/Sm) 



dxx^ ^f{x) 



(30) 
(31) 



with Cf{x) — c^^'^p{cx), using (21 1. Since / is not normal- 
izable, it is not a probability, hence the notation. The normal- 
ization integral 
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of j dSSP{S) converges at small S and 
does not depend~ori the small-scale cutoff |90| However the 
integral j dSP{S) is divergent if we extend the form ( 30 1 to 
small S. It means that 



{S) ~ ^S^i 



(32) 



with fc = (1 — r)/(C/(0)). This is an approximate value, 
obtained setting /J" dSP{S) = 1 and extending the form 



(30 1 down to Suun- Note that although (S) = JdSSP{S) is 
a convergent integral at small S, its actual value depends on 
the precise cutoff at small scale. This is because, for r > 1, 
almost all avalanches are in size of the order of the UV cutoff 
5min- Despite that fact however, all moments (S'p) with p > 
r— 1 are controlled by rare avalanches of size ~ Sm, the large- 
scale cutoff. For r > 1 it is only this part of the distribution, 
i.e. for S ^ S'min, which is universal, up to a multiplicative 
constant. 

Finally, in the case t — 1 both phenomena are present. We 
will not study this case here. 



IV. CONNECTION WITH FRG FUNCTIONS 

The hypothesis of a finite density of shocks implies that the 
functions C'^"-'>{vi, . . . , Vn) and C''-"^(wi, . . . , u„) are contin- 
uous and have no ambiguities at coinciding points; i.e. choos- 
ing a given order for its arguments one can take the limit of 
coinciding arguments, and the result does not depend on the 
chosen order 

The derivatives of C'^"''(ui, . . . , u„) however are distribu- 
tions and they do contain information about shock statistics as 
we now show. 



A. Cusp 

Let us start with the second moment ([12]). Taking deriva- 
tives w.r.t. wi and W2 one finds: 

- m-^L^'^ll"{wi - W2) = u'{w^)u'(w2) + 1 



i 

+ L^'^'^Y S^Sj5{wx - Wi)E{w2 - Wj) + 1 , (33) 



where we have used the decomposition in shocks (I61 and 



separated in the double sum the contributions from the same 
shock and from different shocks. 

Hence the second derivative is the sum of a (5 singularity 
and a smooth part: 

- A"(u;) = m^L-"^ J dSS^piS)S{w) 

+m^L-'' j dSidS2SiS2Pc{Si,S2,w) , (34) 



where 



Pc{Sl,S2,Wi - W2) 



<5(S'i - Si)6{S2 - Sj)S{wi - Wi)S{w2 - Wj) . 

(35) 
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Note that p{S) and the connected joined 2-shock size-density 
Pc{Si, S2) are proportional to the size of the system. Integrat- 
ing around zero yields ||9T| : 



2A'(0+) = m^L-^ j dSS'^piS) 



{S) 



(36) 



using ( [2T] i. This provides a rather nice interpretation of the 
FRG function A{'w) in terms of shocks: the cusp gives di- 
rectly the second moment of the shock size, and the remain- 
ing smooth part of the function gives the correlation of the 
sizes of two (distinct) shocks at different points. A non local 



generalization of this formula is given in (F12 1. 

Let us now verify the scaling for small m. From the above 
one finds: 



2A'(0+) = cm'^S.. 



(37) 



with c c ^ ds ^p{s) — 

/o°°dss2-^/(s)//p°°dssi-^/(s) (the latter for t > 1 
only). This has the expected fixed-point scaling ( 14 1 with 
— A'(0+) — cc/{2Ad)- Below, we will make the convenient 
choice c = 2 to fix the arbitrariness in definition of Sm- 

This correspondence can be extended to higher moments. 
For illustration, let us consider the third moment: 



U'{wi)u'{w2)u'{w3) 

L-^'^ ( dS p{S)5{wx - W2)5{wi - W3) 



+ 2pcrmj 

L^^'^ X! SiSjSk5{wi-Wi)5{w2-Wj)6{w3,-Wk) 

(38) 



From this expression one can integrate 
JJ" dull dw2 /q"" dw3 and obtain 



[u{w) — u{Q)Y 



-3d 



I 



dSS^p{S)w + 0{w^) . (39) 



This can be generaUzed to any order, as discussed in appendix 

Note that in dimension d — the field 'm?{w — u{w)) iden- 
tifies with the velocity field of a 1 -dimensional fluid which 
evolves in time t = according to the Burgers equation 
with random initial conditions. The linear cusp of the third 
moment, Eq. ( [39] l is the analogous for Burgers turbulence to 
the famous exact result of Kolmogorov for Navier Stokes in 
the inertial range, both models exhibiting similarities. We will 
not discuss these connections further here, see Ref. 162] l63l . 
but since the manifold problem is a d > generalization of 
the Burgers equation, we will call Kolmogorov moments mo- 
ments such as ([39]l and their generalizations. 



B. Kolmogorov cumulants and generating function 

We now generalize the results of the previous section and 
construct a very useful generating function which resums all 
shock-size moments. 

Since we are interested in the jump sizes, it is convenient 
to define the Kolmogorov cumulants (by analogy with the fa- 
mous K41 paper |82|, as discussed above): 



[u{w) - w - u{0)]" =m-2»i-("-i)'^i^(«)(w) (40) 
K'-^^w) = 2{A{0) - A{w)) , (41) 

where all if^"^ have a large-L limit, and K'^^^ = 0. We find 
that they are proportional to moments of the shock-size distri- 
bution: 

K^''^w) = m^''L-'^ [ dS'5"p(S')|u|(signw)" + 0(u;2) 

(42) 



(S) 

Note that the leading small-w coefficient, the cusp, is 
the same for [u{w) — w — u(0)]" , [u{w) — w — u{0)]" and 
[u{w) — u(0)]". We thus study the generating function: 

^ \n 

QXL''luiw)-n,-uiO)] - 1 = ^ —L'"^lu{w) ~W~ u(0)]" 

n=2 ^' 

00 , ^ ^ 

= T.^ dSS-p{S)\w\{signwr + 0{w'), (43) 

n=2 ■ 

The following generating function hence has a finite large-i 
limit: 



G{X) := L"''(e^-^''["('^')-«'-«(o)] _ 1) 
= L-'^ J p(S')(cosh(AS') - 

+L-'^ J p(5)(sinh(AS') - XS)w + 0{w^) (44) 

For positive w > this yields: 

G{X) = Z{\)w + 0{w^) (45) 
^(A)-^((c^^)-l-A(5)) (46) 

Hence by computing the cusp of G(A) one has direct access to 
the characteristic function of the shock-size distribution. The 
formula ( |44| is also derived in appendix|A]by another method. 
In the following we reserve the notation Z{X) := A + Z{X) 
and the notation Z for the rescaled version of Z (see below). 



V. TREE AND IMPROVED TREE CALCULATION 

We now compute the cumulants C{wi, . . . , w„) and from 
them, the Kolmogorov cumulants and the generating function 
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G{X). We use several methods which give equivalent results: 
a calculation using replicas, one without replicas and a graph- 
ical representation using static or dynamic graphs. 

In this Section we obtain the form of G{X) to lowest order 
in the e expansion (i.e. zero-th order). It is essentially a tree- 
level calculation, although some loops can be incorporated, as 
we discuss. In a diagrammatic language, one defines suitable 
trees where the building blocks contain loops, which are re- 
summed here. Since we restrict to a uniform w, i.e. to zero 
external momentum, all trees carry zero momentum q = 0. 
This is why the calculations to this order look very similar to 
calculations on a d = toy model. 



To perform the Legendre transform one must invert the re- 
lation ni^Wa/T = r(.j.j,p(?ia), i.e. find the function Ua{w), 
solution of 

Ua{w) = Wa+ R!{Ua{w) - Ub{w)) . (51) 

b 

This solution can be expanded in the number of free replica 
sums. 



A. Method using replicas 



which is also an expansion in powers of R. Inserting in ( 5 1 
generates recursion relations: 



The cumulants C are contained in the generating functional 
Wlw] for connected correlations (for more definitions see Ref. 
1621 ). We focus on a uniform Wx — w, in which case W[w] — 
L'^W{w) and 



W{w) 



in ■r-^ 2 



2!r2 ^ 

ab 



R{Wab) 



OO 



^ p\TP 



E 

ai...ap 



S^P\wa,,...,Wa,), (47) 



where Wab = Wa — Wb- It can be obtained via a Legen- 
dre transform from the effective action per unit volume, i.e 
r(u) = L^'*r[u] for a uniform field Ux = u: 

2 

W{w) + T{u) = ^ XI """^^ 



(48) 



III s—^ 2 



1 



2!r2 ^ 

ab 



R{Uah) 



OO 



1 



p=3 ai...ap 



The F-cumulants are naturally computed in the field theory as 
the sum of all 1 -particle irreducible (IPI) graphs, as in ||36l . 
Each S'^P^ can be computed in an expansion in powers of R 
(at T = it is the usual loop expansion). They have the prop- 
erty that to lowest order S^^^ ~ R^ and S^'p^ ^ RP. Hence 
near d = 4 where i? - e, S'^p) - e^. Note that the S'^p^ 
have a different counting, e.g. S'^'^^ ~ e^. An important prop- 
erty, arising from the Legendre transform, is that W{w) can 
be written as a sum of tree diagrams with vertices made of T. 

The improved tree approximation consists in setting the 
higher cumulants of the effective action T to zero, S'^p^ — 
for p > 2, and then perform the Legendre transform to obtain 
the effective action W{w) and S^p\ i.e. to use 



2!r2 E 

ab 



R{Uab) ■ (50) 



It is "improved" as one keeps the exact two-replica part R, 
which itself has a loop expansion. For instance, from the 
above discussion one has T = Ftroo + 0{e^). 



1 



(53) 



u 



(2) 

a.a\.a2 



Syma„a.^"(^aaJ(^i!l, -<,aj (54) 



,(3) 



1 



-1 



R"'{Waa,){u(^l - 4V,oJ(41 - 4V,03) 



(55) 



This is easy to automatize using Mathematica, the combi- 
natorics being similar to expanding R' (Ylp^p^^)^ W\\h the 
additional difficulty of attributing new labels to repeated in- 



dices. Once the terms in the expansion (52 1 are known, since 
m?Ua/T — W'{wa), one obtains the derivatives of the energy 
cumulants: 



5(p)K,,...zi,J (49) 41,...a„H 



1 



(56) 



(n- 1)! 



J . (57) 



Upon further derivation, we obtain the cumulants of the force, 

C'(")K,...u;„) = {-lYd^,...d^J^-^\wi,...Wn) as 

C'("+i)K,u;oi,...,«ioJ 

= - . . .a„„„4:^„...,„(z«) (58) 

We display the obtained explicit forms for the third and fourth 
cumulants: 



3 
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S^^\wi,W2,w^) = A^sym^^sR'MR'iwis) (59) 



{WI,W2,W3,W4) = ^Symi234^12(-R'l3 " ^23)^14 

(60) 



for the energy, already given in ||62| . We use the shorthand 
notation i?,i2 = R{wi2)- They have a simple graphical repre- 
sentation in terms of tree diagrams with R vertices. For the 
corresponding force cumulants this gives (with everywhere 
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A{u) = -R"{u)y. 



6 
12 



Synii23A'(wi2)A(wi3) (61) 



+Ai2A^3A'i4 + 2Ai2A'i3A;j4 
+Ai3Ai4A'/2 - Ai3A24A'/2 



with Ai2 = A(wi2) and so on. Note that this expression for 



the third force cumulant was tested against numerics near de- 
pinning in Ref |66| (within the improved tree approximation. 



the relation (61 1 is the same for statics and depinning with, 



however, t//j^erenf functions A(w) in each case). 



The present method allows to compute significantly higher 
cumulants, using mathematica. We do not display the full 
form of C''-'"\wi, . . . , w„), which are very tedious, but we 
give the resulting Kolmogorov cumulants: 



^2(A(0)- AH) 
m^K^^^w) ^ -12A'(u>)(A(0) - A(w)) 

m'^K'-^^w) = 120(A(0) - A{w))A'{wf - 48(A(0) - A{w)fA"{w) 

m'^k^^Hw) = -80(A(0) - A{w)){2lA'{w)^ - 24(A(0) - A{w))A'{w)A"{w) + 2(A(0) - 

m^k'^^^w) = -480(A(0) - A{w)){-63A'{w)'^ + 138(A(0) - A{w))A'{wfA"{w) 

-22(A(0) - A{w))^A'{w)A"'{w) + (A(0) - A{w)f{-18A"{wf + (A(0) 



A{w)fA"'{w)) 
- AH)A""H) 



(62) 
(63) 
(64) 
(65) 

(66) 



We expect, from the assumption of a shock density (see the 
discussion above), that the C*^"' are continuous functions of 
their arguments. This we checked explicitly. Hence there are 
no ambiguities at coinciding points, i.e. to perform the limits 
one can choose any order for the arguments, with the result 
being independent of the chosen order 

One first checks that the values at zero vanish: 



{u{w) - w)" = u(0)" cx (0, 0, . . . , 0) = (67) 

Hence the distribution of the center-of-mass position in the 
quadratic well is Gaussian to lowest order in e, i.e to the im- 
proved tree approximation. To this order the calculation is the 
same in statics and dynamics. As we see below corrections to 
the Gaussian arise to one loop order (see |51 , 83| for some re- 
sults on the deviations of the distribution of the interface width 



to the Gaussian at depinning). At quasi-static depinning (67 1 



also gives the cumulants of the distribution of the critical force 
ll64l . Hence it is also Gaussian to this order This is consistent 
with Ref. 153 1 were deviations from Gaussian were found and 



computed to one loop. In fact (67 1 validates these calculations 
within the present well-controlled setting of a quadratic well. 

We now expand the above result for k^'^^w) to small ar- 
gument and find: 



k^^\w) = a„(-A'(0+))"-ii«(signu;)"m^ 



4-2n 



(68) 



with 



02 = 2 , 03 = 12 , 04 = 120 , 05:= 1680 , 

06 = 30240,... (69) 

Note that o„+i/a„ = 4n — 2. Defining &„ := a„/n!, the first 
coefficients are 

62 = 1 , 63 = 2 , 64 = 5 , 65 = 14 , 66 = 42 . (70) 



B. Method without replicas 

An equivalent method to this order is as follows. One no- 
tices that each cumulant C*'"' is computed to lowest non- 
vanishing order in R (or A). Hence it is formally equivalent to 
start from a replicated action S containing R only (i.e. a bare 
disorder with the substitution Rq — > R) and compute the mo- 
ments of u{w) — w each to lowest order in perturbation theory. 
Hence the calculation of the previous paragraph is equivalent 
to the following one in dimension d = 0: Denote u{w) the 
minimum of the toy model: 

^2Htoy=liu^wr+^^V{u), (71) 

where V{u) is a Gaussian random potential of correlator 
R{u). The minimum satisfies: 



u{w) — w -\ -F{u{'w)) 



(72) 



with F{u) = —V'{u) and A has been introduced to count the 
powers in V. Compute each moment defined as: 

(-l)"m-2"i-("-l)d(7(")(u;4, ...,Wn) 

{u{wi) ~ Wi) . .. {u{Wn) - Wnf 

^ ^ F{u{wi))...F{u{w„)f (73) 

perturbatively in A (in fact in \/m?) to lowest non-trivial or- 
der, which is 0(A^"^^). The factor ^^-("-i)'' can be omitted 



since d ~ Q. One uses iteratively ( 72 1 and the Wick theorem 
with F{w)F{w') = A{w — w') with all other cumulants of 
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F{'w) set to zero (i.e. a Gaussian F). Using Mathematica we 
have reproduced most of the resuhs of the previous Section. 

Note that this is different from the standard perturbative 
expansion of the toy model which yields dimensional re- 
duction, i.e. a trivial perturbation expansion involving only 
A(0) — — i?"(0). The difference is that one computes here 
cumulants (7*^"'' at different points and only at the end the limit 
of coinciding arguments is taken, using a non-analytic A{u) 
with a cusp. Nevertheless ambiguous terms at intermediate 
stages may be generated. To the lowest order studied here, 
those represent possible contributions to C''^"-' which depend 
on a smaller number of points and they cancel in the calcula- 
tion of the connected correlations (7^"^ 



C. Graphical representation 

It is useful for the following to give a graphical representa- 
tion of the results of the two previous sections. Define: 



Wi, 



,Wn) 



For instance we write the 2-point correlation as 

1 



i-A(wi - W2) 



(74) 



(75) 



The graphical notation here and in all diagrams of this type 
are as follows: there are external legs with points on the top 
labeled by integers 1 to n, corresponding to positions w\ to 
Wn and external fields u(wi) ~ Wi, here n — 2. The A 
(or equivalently R) vertices are double vertices (non-local in 
w) with two points joined by a dotted line, and can be inter- 
preted equivalently as in the statics or in the dynamics (we 
have checked equivalence to the order we are working). In the 
statics, they are R vertices, giving ~R" after contraction with 
the u-fields. In the dynamic formulation, they are A vertices 
and the two lines exiting a vertex are directed to the top, and 
end up being equal to a static propagator: in the real dynam- 
ics they are response function, which usually are denoted with 
an arrow: we do not show the arrow here but they are always 
implicitly towards the top of the picture. Here they are always 
evaluated at zero frequency. Thus the lines are static propa- 
gators, evaluated here all at g = 0, hence giving a factor of 
1 /m^. For a generalization to non-zero external momenta see 
Appendix |F| 

With these diagrammatic rules there is a single diagram to 
represent the third cumulant: 



C(3) 



6 



Sym Ai2A'i 



3 ' 



(76) 



where the combinatoric factor comes from the 6 inequivalent 
ways to assign three labels to external legs. The result agrees 



mulant: 

C(4) = i- 





1 




1 


\ 1 


» 1 


L::::: 




I 





(77) 



which reproduce each term in (61 1, all with a factor 4! = 24 
and a factor 1/2 from the three diagrams symmetric in ex- 
change of a pair among the four labels. 

To compute the n-th cumulant of {u{w) — w — u{0)), 
we must evaluate C*-"^ for Wi w minus Wi 0, 
for each i. Writing this operation as ICf{wi, . . . ,Wn) '■— 
riLi ^w)- {w., -> 0)] f{wi, Wn), this gives 

if(2)(w)^/c[ [ =2[A(0)-A(w)]«-2A'(0+)|u;| 

(78) 

where from now on we extract the m dependence from the 
lines of the graphs, which are hence set to one. Note that 
there are four choices to assign or w to each leg, hence four 
terms. For the third cumulant one finds: 



-1 



-12[A{0)- A{w)]A\w) 
l2A'{0+fw 



(79) 



which produces 6 x 4 = 24 terms. 
Similarly one finds: 

^ 120 [A(0)- AH]A'(w)2 
« 120A'(0+) \w\ 



48 [A{0)-A{w)fA"{w) 
(80) 



in terms of the five diagrams in (77i. These results agree 
with those of the two previous sections. Here too the con- 
tractions are unambigous, as long as all Wi are different. It 
would appear naively that it is equivalent to apply JC to an un- 
symmetrized expression or to its symmetrized form. This is 
not true in fact because of ambiguities at coinciding points. 
One must be very careful to apply K, to the symmetrized ex- 
pression and then take the limit, i.e. 



(81) 



lim 



with (61 1. Similarly there are five diagrams for the fourth cu- 



n [{w^ ^w + S,)-{w,^ S'^] C^") . . . , w„) 

i=l 

where the Si , 6'^ are taken to zero in such a manner that all 
arguments remain distinct and the order is fixed, the C*-"^ be- 
ing fully symmetric functions of their arguments. We have 
checked here and in the 1-loop extension given below that the 
result does not depend on the order as required by the conti- 
nuity of the force correlations. 
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D. Recursion relation and resummation 



It is clear from the above results (|78]l-([80j and ( |66) that 
the leading term in the expansion of K^^\w) in small w is 
obtained from those terms which contain a single factor of 
A(u)) — A(0), since the latter is of order w. There are also 
terms such as [A(i(;) — A(0)]p with p > 2 in the expressions 
for but they are of higher order (since e.g. A"(w) 

A"(0+) has a finite limits at w) = 0+). Therefore, the final 
result ( [68] l for the coefficient of the cusp depends only on the 
part of C'"' of the form [A(u;) - A(0)] A'(u;)""2, and is the 
diagram with one and only one terminating (lower) A vertex. 
This is not the case e.g. for the third and fifth diagrams in (j77]l 
which have p ~ 2. 

Let us call C'"-* := (— l)"^^") the part which contains only 
one factor of A{wi) for some i (we do not write explicitly the 
— A(0), but remind that every A{w) comes with a — A(0)). If 
we restrict to this part we can write a recursion relation using 
the tree structure. Indeed, to construct C^"-* we can either: 

(i) takeC(""i) and a single (as yet unlabeled) leg u{w) — w, 
and attach a A vertex at the bottom (hence without deriva- 
tives) to each of these two elements. It results in a diagram 
with again only one lower A, e.g. 



-i 



-I 



(82) 



and when adding the two new terms I = 1 and I = n — I, the 
combinatorics is identical. We can indeed take advantage of 
this feature. To do so, we go to the Kolmogorov-cumulants 
and set as in (68i and for w > 0, m^"~^i^'"'(w) ~ 
a„[-A'(0+)]"-2[^(Q) _ m2"-4ii:(0'(o+) = 

ai[-A'{0+)Y-\ We note that gluing X(') and yields 
i:(')'(0+)ii:("-')'(0+)[A(0) - A{w)] times the combinatoric 
factor ( " ) . Indeed one glued part should be at Wi = w and 
the other one at w = 0. Hence we can convert ( (84]i in to a 
recursion relation for the a„. Inserting the ansatz (|68]l, the 
above recursion relation becomes: 



ji-i 



1=1 



°-n^^{ V ) aittn^i , ai = l . (85) 



We note that it correctly reproduces the series (69 1. The series 
bn ■— a„/n! then satisfies the recursion relation 



(86) 



Hence the generating function 



Zo(A) := 



(87) 



The combinatorial factor is n (here n = 4) for labeling the 
newly added leg. 

(ii) starting from the fourth cumulant, take two already con- 
structed cumulants (trees) of size I > 2 and {n ~ I) > 2, and 
glue them together using again a A vertex at the bottom, e.g.: 



(83) 



Now the combinatorial factor is ( " ) , for choosing the group 
of indices for the first and second element, with the restric- 
tion that I < n — I. If n = 1/2, an additional factor of 
1/2 appears. It is more systematic to sum over all pairs, i.e. 
I = 2, ... ,n — 2, and divide this sum by 2. This can be 
summarized in the following schematic recursion relation for 
C^"^(wi, . . . ,Wn), which by construction is a function of the 
set of Aij — A{wi — Wj) and A'^j = A'{wi — wj): 



f-{n) 



P)r(n-1) 

AUA 



i<j 
n-2 



<5A,: 



) n 



i=2 i<2 ■' k<m 



ij km 



X (A,,fe + A 



3,m 



We have temporarily set the mass 77i = 1 to not burden no- 
tations. We remark that the first term has formally the same 
form as the others, except that one of the derivative-terms does 

not exist. Since n = ( " ) 





1 








~ 2 




-(„=.)] 



satisfies the equation: 

Z{\) = X+Z{\f . 

It has solution: 

1 



Z(A) = Zo(A) :=-(l-Vr^) 



(88) 



(89) 



Let us now recall that from ( |40| i and ( |68| ) one has at small w: 
1 



n 



{w)-w- u(0)]" 



(90) 



A'(0+) 



ri-l 



w;(signH)"-^ + OK) 



Following the discussion of Section IV B|we find that: 

(i) the generating function defined in (|44[), ( |45] l for w > 
from G(A) = Z{X)w + 0{w^), satisfies the scaling form 



Z{X) = -^Z{XS^) - X 



where from now on we define (recall A'(0+) < 0) 



|A'(0+)| 

m 4 



(91) 



(92) 



Since (36i is exact, Eqs. (j9TJ and ( [93] l are valid to all orders, 
as they amount to the choice Z{X) = A + A^ + 0( A'^), i.e they 
fix the coefficient A^ in Z to one. Scaling means that in the 
limit TO — > the function Z becomes m independent, whose 
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validity is discussed below. In summary, using the exact rela- 
tion ([36]l the definition of Sm chosen in this paper is: 



(93) 



(ii) in the (improved) tree level approximation the function 
Z{X) satisfies the self-consistent equation (881, with solution 



Zq{\) given in (89 1, but with A'(0) renormalized. Hence we 
find the generating function of the avalanche size distribution: 



:-7^[(c^^)-l-A(5)] 



1 

Vs) 
1 



(94) 



with S„i given by ( |93| . Note that the presence of the factor 
1 / (S) indicates that one can only obtain information about 
the distribution P{S) up to an unknown multiplicative factor 
(S), consistent with the discussion in Section (lllC 1. We now 
analyze this result. 



E. Moments of avalanche sizes and universal ratios 



The easiest quantities to extract from ( |90| ) are the moments 
of the size distribution P{S). From (94 1 one finds for n > 2: 



(S) 



n-1 



(95) 



The coefficients a„ and 5„ can be calculated explicitly for n > 
1: 



d_ 
dX 



1 



1 - Vl - 4A ) - A 



A=0 



(2n- 3)! 



(96) 



n! (n-2)!n! 
Although there are some cases where the scale Sm given by 



(93 I can be calculated, see below, usually it contains a non- 
universal amplitude. Hence it is interesting to form universal 
ratios independent of any scale, such as: 



(97) 



One finds their value at tree level: 



F. Distribution P{S) 

Let us now perform the inverse Laplace tranform of Z{X) to 
obtain P{S). One must be careful since one expects a form for 
P{S) which is not normalizable in the absence of a small-S" 
cutoff, hence we work with SP{S) / (S) which is normalized. 
From ([94| one can write (setting momentarily Sm ^1) 



(S)-^ / dSP{S){e 



AS 



1) = -(1 - ^/^^4A) , (100) 



where the S integral converges at small S. This is equivalent 
to 



{sy 



dSSP{S)e 



AS 



1 



VI -4A 



(101) 



as one can check by integrating over A on both sides. Using 



'too 27r? 

{S)-^SP{S) 



5{s — sq) one has by inverse Laplace 



dA 



-AS 



^oo 27ri Vl - 4A 

TT 7l/4 V4A-1 



g-AS^ 



(102) 



^-1/2 



upon deforming the contour around the branch cut. Restoring 
all factors yields the final result: 



PiS) 



{SI 
2^ 



^-l/2^-3/2g-S/(4S„ 



(103) 



As discussed in Section ( 111 C ) this is expected to be valid for 
S ^ S'min in the limit of small m, large Sm- Note that in 
general the exponent r can be extracted from the behaviour of 
Z(A) for A —00, which is dominated by small avalanches, 
and reads, in the scaling regime: 



zix) = x + zix)^ -s;,-'\x\ 



(104) 



On the other hand, from its definition (|45| and ( |9T] ) we see that 
A + Z{X) must converge to -1/(5) - -S^T^^S';-^ at A ^ 
—00; thus the crossover occurs for A ^ 1 /S'min- This estimate 



is consistent with the relation (32i. For larger values of —A, 
Z{X) becomes non-universal with respect to UV details and is 
out of reach of the present method. 

We can now compute the moments for arbitrary real n > 
1 /2 by direct integration from ( 103 1. One finds: 



i2ri-2 



r n 



(105) 



which agrees with the analytic continuations of the moments 
obtained above. 



f 2 = 3 , 



2n-3 ' 



T'3 



(98) 



(99) 



G. Discussion of tlie result: mean-field theory 



The results (|93]l and ( |103| for P{S) imply the value of the 

(106) 



avalanche distribution exponent 

r = 3/2 
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FIG. 1 : Example of a diagram at MF level, as generated by Eq. ^142\ 
at a = 0. It contains a correction to disorder i.e. A'(O^) at 1 loop 
(shaded in gray). 



It agrees with the results obtained in: (i) a mean-field toy 
model for dynamic avalanches at depinning fSl, (ii) a mean- 
field argument given in the context of a non-equilibrium 
random-field Ising model 122, (iii) mean-field calculations 
developped for sandpile models llT4l [15 J . This form hence 
seems rather robust as a mean-field result. 

Here however, it is derived from first principles using FRG 
from the elastic model ([TJ. It was obtained by resumming 
all "improved tree" diagrams, i.e. trees made of fully dressed 
vertices R. In the standard bare perturbation theory it includes 
diagrams as represented in Figure[T] Hence it can also predict 
the scale S„i given by ( [93| ). 

For the elastic model, the result ( [93| l, ( 103 i for P{S) is 
valid for dimensions d > 4. Convergence to this result as 



m ^ requires a 1-loop analysis, as discussed below. The 
only point to discuss is thus the scale 5„i, i.e. the m depen- 
dence of A' (0+) . For that we use some of the discussion about 
elastic manifolds for o? > 4 and the associated (Wilson) FRG 
flow, given in Appendix H of |42|, and some specific results 
for the present model, summarized in Appendix[B] The rough- 
ness exponent is = 0, and the manifold has a finite width. 
If bare disorder is smooth, it should be sufficiently strong for 
metastability and (typical) shocks to exist (a cusp in A{u) then 
develops). Alternatively one can consider weak rough disor- 
der (i.e. with a cusp), or smooth but with a very short correla- 
tion length. In these cases one has large avalanches with 



|A*'(0+)| 



(107) 



where for d > 4, |A*'(0+)| is a non-universal number, see 
Appendix [B| One can compare this scale with the fluctuations 
of the center of mass, such that L'^u{0)^ = A*(0)/m^, also 
non-universal. In c? = 4 one recovers universality and finds: 



87r2|A*'(0+)|m-* 



In 



Too 



i-Ci 



(108) 



at small to, while, for comparison, L u{0)^ = 
A*(0+)m-4[ln(^)]i"2Ci^ where Ci and aU fixed-point 
values are given in appendix [B| 

In the end it may not be surprising that mean-field theory 
is obtainable as a summation of trees. However, this had not 
been done previously within the FRG. It is quite remarkable 
that the result for the mean-field generating function Zq{X) in 



rooted binary planar trees with n + 1 leaves 1*921, also known 
as the Catalan numbers C„ = 6„+i. The latter also appear for 
the number of rainbow-diagrams in RNA folding, the number 
of Dyck words of length 2n, the number n + 1 factors can 
be completely parenthesized, the number of monotonic paths 
along the edges of a grid with n x n square cells, and many 
more. It would be interesting to obtain a more microscopic 
understanding of the relation between the diagrammatic trees 
and the structure of the (static or dynamic) avalanche pro- 
cesses in d > 4. 

This task may be easier to carry out for the dynamics. Intu- 
itively an avalanche starts with a seed and each event may or 
may not trigger new events. If the latter do not influence much 
each other, as expected in mean field, the process indeed looks 
like a tree. These features are captured by the simple model 
studied in Ref. [8|. There, an avalanche starting at time < = 
and lasting until t = T, of total size 5* = X^t^o "^t ^ succes- 
sion of jumps of size 5'^ each being the sum of rij independent 
events, i.e. St = ^i- The number of events rit has a 

Poisson distribution of average {fit) = pSt-i, i.e. depends 
only on the jump at time t — 1 (with tiq — 1, and St — for 
t > T). If the sizes of the independent avalanches are simply 
Si — \ this model is the famous Galton process |84|, which 
describes the evolution of a population of size St with pois- 
sonian distribution of numbers of offsprings. In this model 
criticality arises as p is increased up to the threshold for an 
infinite avalanche. Below the threshold the distribution of the 
avalanche size S is exactly the one obtained here by summing 
the trees. The S~^/'^ power law distribution may then be un- 
derstood as the distribution of return time to the origin of a 
random walk. Interestingly, there is also a d = model which 
captures this mean field physics, and its connection to return 
time of a random walk. It is the celebrated ABBM model ||85l 
for domain wall motion, represented as a particle in a Brow- 
nian random force landscape. We recently computed 1861 the 
avalanche statistics and the renormalized FRG force correla- 
tors for this model, and the similarities with the mean field 
results obtained here are striking. These analogies involve 
mostly the zero momentum, q = Q structure, but in Appendix 
|F] we present a spatial generalization of our tree summation, 
which may help to understand the relation between trees and 
avalanches in mean-field. 



VI. LOOP CORRECTIONS 



A. General method 



Until now the force correlations (7^"^ and from them the 
K^"-^ and the avalanche size moments {S^)/{S), were com- 
puted to improved tree level, i.e. setting S*'"^ = for n > 3 in 



( 89 1 is identical to the generating function of the number of 



the effective action ( 49 1. After Legendre transform this gives 
each C**-"' to the lowest order in an expansion in powers of A 
(or H), and later, using the fixed-point values for A, to lowest 
order in e. 

It is useful to describe a systematic procedure to compute 
these quantities to higher order in the expansion in powers of 
R in the statics (or A in the dynamics). For each specific 
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calculation, an intuitive diagrammatic representation is given 
below. 

(i) one first computes all the functions 5"^"^ n > 3, in the 
effective action in an expansion in powers of R |93|. They 
are given by the sum of all 1 -particle irreducible (n-replica) 
diagrams with R as vertices, i.e. the sum of all loops. For 
instance, to leading order to which we restrict below one finds 
S'l") ^ i?" with a unique 1-loop integral /„: 



(fc2 + m2)" 
1 

k (fc2 + l)" 



(109) 
(110) 



for d < 2n, and only the replica combinatorics is non-trivial, 
e.g. for n = 3,4: 



sill = 3/3sym[(R'/2)2R'/3 - ^R'{,R'^,R'^,] (111) 



1 , 



S[%, = 3/4sym[2(R'/2)2R'/3(2R'/4 + R'^^) 



A D" D" D" D" 

-4K12K13K23K14 



R" P" P" P" 1 



(112) 



where here R"{w) := R"{w) - R"{0), R'^^ = R"{wa - 

Wb) and s[^T^ = S''^^^(wi, u'2, 103) and so on. The formula 
for general n was given in I35J and takes the form for the n- 
replica term: 



c(r 



{n-iy. 



In 



tr(iy")-5]i?L^...i?: 



■/f 

'aa„ 



(113) 



i.e. it is a trace over replica indices, with the n + 1 replica term 
subtracted. This formula yields \\\\\ and ( 112i, the explicit 
form for n = 5 was given in ll42l . 

(ii) one then performs the Legendre transform ( 49 1 
from S-^") to S-f"). From them one can then obtain 
C(")K,...,«;„) = ...9^„^(")K,...,?i;„). 
Legendre transformation simply means that each S**^"^ is a sum 
of all (n-replica) tree diagrams which can be drawn using R 
or any of the 5*^"^ n > 3, at the vertices of the tree. Since 
when forming a tree the number of replica can only increase, 
there are simple exact formulae for the lowest moments: 

^^'H«^123) ^Sf\wi2z) + 5(^)^23) (114) 
(U'1234) + '5'('^)(wi234) 

+ ^symR' {wu)d^,S^^\wi2^) , (115) 



as can be checked by explicit Legendre transform 11621 WS 



Here the 5, 



(") 



i?" ^ are the result of the improved tree 



approximation, i.e. the trees made of only R at the vertices. 



as described above, giving (59 1 for the two lowest moments. 
Note that S'^'^) contains a single loop integral while 5'-'*) 
contains both 13 and I4, and so on. Since everywhere we 
consider uniform w, the trees are at zero external momenta 
(hence the factor 1 /ni^), and loop momenta only flow inside 



each S'(") vertex. In principle using ( 115 1 and inserting each 
S'") to the required order in R allows to compute the C^") to 
any desired order in A and in e. 

Here we only want S'(") to order i?", hence we do not need 
all trees made of S'^"^ and R, but only those trees containing 
one S'(P) and one or more R, i.e. it with schematically the 
form: 



in) 



n— 3 



sym. 



(116) 



We first perform an explicit calculation of third and fourth mo- 
ments to one loop, before attempting the resummation. Note 
that to 1-loop order, an equivalent procedure is to perform per- 
turbation theory in the bare disorder Aq, in which case an ad- 
ditional term I2 will appear in each C'- '^K with trees based on 
a single Aq loop. Inserting the 1-loop expression of Aq as a 
function of A this term disappears and one recovers the same 
result as with the above method. 

We need below the momentum integrals /„. One has /„ ^ 
ds exp(— s — sk'^)s"^^ /{n ~ 1)!, hence: 

^ T{n - d/2) 
e/2 2r(n)r(3 - d/2) 

rf^4 r(n - 2) _ 1 



2r(n) 2(n-l)(n-2) 



(117) 



where e/2 is discussed in Appendix |B| Note that all these mo- 
mentum integrals are IR-finite for d < 6. This will remain 
true to any number of loop as all diagrams entering in C**-") 
are superficially IR convergent by power counting for d < 4. 
If one uses perturbation theory in the bare disorder, they will 
however contain diverging subdiagrams, starting at 2-loop or- 
der These divergences should be removed by the countert- 
erms for the disorder, i.e. the replacement of Aq as a function 
of A. Therefore, corrections to r should come from finite dia- 
grams, and appear when summing all diagrams. On the other 
hand, for small n, the moments are interesting since they can 
be compared to the numerics. 



B. Third cumulant 



From the above expressions ( |lll[ i and ( 1 15 1, we obtain 



6hsym[A'{wi2fA'{wi3) + {A{wi2) - A(0)) 

X {A'{w,,)A"{w,2) + A'iw,2)A"{w,,) 
+ A'iw23)A"{wu))] (118) 

It is useful to indicate the graphical representation of each 
term, using the conventions described in Section |V C| (sup- 
pressing for clarity below all propagator and momentum loop 
factors). There are three non vanishing 1-loop contributions: 



-6 Sym A'/2A^3 (A12 - A13) h 



(119) 
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= 6(Ai2-A22)A'i2A^'3/3 

(120) 



-6[a;2]^a'23J; 



(121) 



They reproduce the above result ( 1 1 8 i after symmetrization. 
One additional diagram 



Sym ;' 



Sym A;2A^3A'3i = (122) 



vanishes after symmetrization in the statics. 

Note that the total result for C'"^) should be continuous and 
unambiguous. Note that C^^) (0,0,0) = while this IS not 
the case for the dynamics (it gives the third cumulant of the 
critical force). 

From there we compute the Kolmogorov cumulant 
K'-^^w). One finds for each piece: 



/C 



12/3 [A(0) - A{w)] A'{w) [A"{w) + A"(0)] 
-24/3A'(0+)2A"(0)u; (123) 



/C 



-1 



-I 



= 12/3(A(u;) - A(0))A'(u;) (A"(0) - A"{w)) 
= O {w') 



(124) 



and 



= 12/3A'(u;) (A'(0+)2 - A'{wf) 
= ^2AhA'iO+fA"{0)w + O (w^) 

Together, they give 



(125) 



K 



(3) 

1 — loop 



-i 



-i 



+ 



Hence, using ( 42 1 we obtain the third moment of avalanche 
sizes: 

(S^) ^ i|^.(o+)2 h - 4m2/3A"(0+)] + 0(A4) 



(5) mS 
= 125^(1 -4A"(0+)4) 



(129) 



We now use that 

e/2 4 ' 



A"(0+^ 



(130) 



where the second equality is valid at the fixed point 
(see Appendix [b|). Using the exact relation {S'^)/{S) = 
— 2A'(0+)/m'' = 2Sm, we thus obtain the universal ratio to 
one loop 



(2) 

1 — loop 



js) js') 



1 - A"(0) 



(131) 



C. Fourth and higher cumulants 



Similar calculations, either directly from ( 1 13 1, or by cal- 
culating the sum of all diagrams at 1-loop order, we find the 
corrections for the fourth Kolmogorov cumulant 
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QA'{wf (A'(0+)2 + 4(A(0) - A{w))A"{w)) 



A'{Q+Y + 2(A(0) - AH)^A"H^ + 7A'{wy 



h 



+ 48 (l6(A(w) - A(0))A"(w) - 3A'(0+)2) A'{wf 

+ 4(A(0) - AH)2A(3)(ii,)A'(u;) + iA'{wf 
+ 2(A(0)- A(w))A"H 

X (A'(0+)2 + 2(A(0) - A{w))A"{w)) ] h 
In the limit of small w. 



K 



(4) 



1-ioopH = (480/4 + 960^ ) A'(0+)W(0+)|z«| 



(132) 



= 12/3A'H [A'(0+)2 - A'{wf - 2{A{w) - A(0))A"(u;)] 
= -m^iA' {Q+f A" {Q)w + O {w^) (126) 

The final result for the small-w behaviour of K'^^^ (w) is: 



,(3) 



-loop 



12 



A'{0+f [l-4m2/3A"(0+)] 



w - 



(127) 
(128) 



Adding the tree-level part, and using ( 42 1 we obtain the fourth 
moment of avalanche sizes: 

^ - ^ l^'(0^)|' [1 - 4 {m% + 2m^/3) A"(0+)] 
= l20Sl 



1-4 



ei2 



(133) 
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It yields for the universal ratio: 

(3) _ jS') jS') 



1 — loop 



(53 



1 - 4^A"(0) 
eh 



1 - A"(0) 1 - 



6 



1 - 



{e-C){6-d) 
18 



(134) 



With Mathematica we have computed C**^"^ and if up to 
n = 6. The expressions are tedious. Let us just note that one 
finds the limits at zero argument: 



C(2p)(0,...,0) = -(2p-l)!/2pA'(0+) 



2p 



(135) 



for even n > A and zero for odd cumulants. This allows to 
compute the probability distribution of the center of mass fluc- 
tuations. The calculation and results are given in Appendix [G] 
It bears some similarities with the one for the distribution of 
the critical force at depinning in ||53]| . 

Let us give here only the part of the Kolmogorov cumulants 
proportional to /„. They read: 

k^^){v) = -240/5A'(w)(-A'(0)4 + 3A'(w)4 (136) 
+ A^{vf^"{vf - 2A'(w)2(A'(0+)2 - 8A(w)A"(w))) 
k^'^Xv) = 240/6(-A'(0)* + 3lA'(w)6 + AA{vfA"{vf 

- 15A'(v)4(A'(0+)2 - 16A(?;)A"(t;))) 

- 15A'(w)2(A'(0+)4 - 8A(w)2A"(i;)2)) (137) 

Expanding at small w > one finds, for the part proportional 
to /„: 



-A„/„(-A'(0+))"-iA"(0)«; + 0{w^) 

(138) 



with coefficients: 



can absorb all factors of m and A'(0"J^by using the rescaled 



generating function Z{\) defined in (91 



and (93 1. One has 



(141) 



where the c„ are extracted from the Kolmogorov cumulants 
at small w, and at tree level c„ = a„. In the previous sec- 
tion we have obtained the contribution to c„ proportional to 
/„, involving the coefficients A„. Remarkably, the remain- 
ing corrections (i.e. the terms in c„ proportional to Ip with 
?> < p < n) can be generated automatically by taking ad- 
vantage of the tree structure. This is illustrated graphically as 
follows: 








Here, the shaded blob represents Z{\). It must contain all sets 
of trees made with R vertices joined together at their basis 
either by a i? vertex or by a S'*^"-' vertex with n > 3, i.e. a 
loop, in the spirit of \\ \6). Expressed as an equation it reads 



Z{X) = A + Z{Xf + a V(n + 1)2'^-^ inZ{\y' , (142) 



n>3 



where we have defined z„ :— /„/(e/2), and substituted the 
factor (n + 1)2"^^ = An/nl. Here and below we denote 



a := -€i2m-'A"{0+) = -A"(0+) 



(143) 



which is of order e. It was not strictly necessary to make equa- 
tion ( |142| l self-consistent, one could with the same accuracy 
of 0(e) replace all terns Z" with n > 3 by . By mak- 
ing it self-consistent we include more loops. It is however 
more convenient in view of the form at tree level. It also has a 
nice interpretation in terms of the generating function of trees 
with higher-n branching: these occur at a small rate, which is 
0(a), and n-dependent. 



^2 = 6 , A3 = 4 
A5 = 5760 , Aq 



, Ai = 480 
80640 , . . . . 



(139) 



Note that An+i/An — 2n + 4. This suggests that 

A„ = 2"-2(,j + i)! (140) 

which is supported by a heuristic argument in appendix [C] 



VII. RESUMMATION OF ALL DIAGRAMS WITH A 
SINGLE LOOP: SELF-CONSISTENT EQUATION 

A. Self-consistent equation 

We can now perform a resummation of all cumulants to 1- 
loop order In view of the results of the previous sections we 



B. Summation of the series and final equation 

From the above self-consistent equation one can compute 
iteratively Z{X) with the result: 



Z{X) ^ A + A2 + 



120 



X^ + 



120 



A3 + 



230400^^^ + 20160m3 



1 r 
24 



960m3 + 480m4 



14400ai4 + 5760m5 + 1680 



A^ 



1 

720 



1290240i^a2 + 806400^3*4^^ + 483840^3^ 



+ 403200i4a + 241920^50 + 80640^60; + 30240 
O (X') 



A*^ 
(144) 
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One recognizes that the results from the third and fourth cu- 
mulants of the previous sections are correctly reproduced. 
Note that the terms with higher powers of i„ contain more 
than one (non-overlapping) loops and higher powers of a = 
0(e). The terms of order a are all correctly generated. 

We now want to compute the infinite sum in Eq. ( |142[ ). For 
this we use the following representation of the i„, as in the 
derivation of ( |117| i: 



with TV = l/(2r(3- I)). One finds: 



(145) 



(146) 



n>3 



^0 i^.Hn 



1)' 



that from \ \A2\ it can also be written as: 

Z{\) =: Z = \+Z'^ 

a_ r r z^ 



z 



z 



fc2 + 1 



- 3- 



(fc2 
Z2 



1 - 2Z) 

(151) 



■(fc2 + l)2 

As shown in Appendix [C] this formula has a simple graphical 
interpretation. The last two terms can be thought of as coun- 
terterms which fix the coefficients of Z and Z^ according to 
the choice we implemented here. (They would be absent if we 
wrote down an expansion in the bare disorder.) 



VIII. 1-LOOP RESULTS FOR THE AVALANCHE-SIZE 
DISTRIBUTION 



A. Avalanche exponent r 



Let us define: 



ri>3 

One can write: 

^(n+l)2"-\,z" 

n>3 

^JV / dse-''s-'*/2-V(s^) 



(147) 



z (^(2 - {d + 2)z)(l - 2z) V + 3dz - 6z - 2^ 



1 

■ 2 
2z' 



{d-2){d-A) 
z(2z + (l-3z)log(l-2z)) 

+ — + 2z-' 



(148) 



The last line is also obtained by using j„ 
Eq. ( 
11^2 

Z{X) 



— — ^ from 

Eq. (|1 17[). One can thus rewrite the self-consistent equation 



X + Z^ 



(2 - (<i+2)Z)(l - 2Z)^ +UZ -&Z -2 
(d- 2)(d-4) 



(149) 



Close to d = 4 it becomes 



1 ~, 



Z=\ + Z^ + a-Z[2Z + (1 - iZ) log(l - 2Z)] . (150) 

This equation has to be inverted in order to get both the expo- 
nent r and the tail of the distribution P{S) for large S. Note 



The formula \150\ contains all we need in order to retrieve 
the avalanche size distribution. We recall that 



5-/^(A5„,) = (^)-^((c^">-l) 



(152) 



from ( [9T] l, where Sm is given by ( [93| . Note however that 
( 150 1 contains only information about sizes of order Sm, from 
5 <C Sm to S* ^ Sm, but in all cases much larger than the mi- 
croscopic cutoff, i.e. S ^ Smin, which has not been included 
(and was not needed) in the above analysis. 

Let us start with the exponent r. It can be extracted as 
follows. There is a critical value pc < 1 such that all mo- 
ments / dS P{S)SP diverge for p < pc- Of course these do 
not strictly diverge, as they should be in the end cutoff by 
Smin, but as just discussed we can forget that here to extract 
the T exponent. This exponent is defined from the behavior 
P{S) ^ S^'^ for Smin ^ S ^ Sm- Hcncc we can identify: 



T ^ 1+Pc 



(153) 



We start from the identity 



n-p) Jo 

jS) 1 

Sm r(-_p) Jq 



dXX-'-P Z{-Sm\) , (154) 



where in the second line we have used ( 152 1. Since Z ^ X at 



small A the integral converges at the lower bound. However, 
it may diverge at the upper bound. Hence pc is obtained from 
the power-law tail of Z{X), for A —oo. The boundary is 
for Z{—X) ^ — AP'=, where we have indicated the sign of Z. 
(Z is negative, due to the dominance of the term —1 in the 
definition of Z for large negative A). This is consistent with 



the asymptotic behavior ( 104 1 and the discussion there. 



From (|150| we obtain, for large negative Z = Z{X) 



Z2 



l + a 



-alog(-2Z) 



i-z) 



2-4q 



0(a2) 
(155) 
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This gives ~Z « (-A)i/(2-f ") ^ (-A)p=, from which we 
identify 



1 3 3 

s — ~ — I — a 

2 -la 2 8 



(156) 



(157) 



Using relation (B12i of Appendix [B] valid for all classes of 
disorder, we obtain our main result for the avalanche expo- 
nent. 



c 



+ 0(6^) 



(158) 



Let us now discuss the significance of this result. First we note 
that this formula agrees to first order in e with the conjecture 



Tc 2 - 



(159) 



As mentioned in the introduction, this conjecture was put for- 
ward in the study of interface depinning in Ref. [271 in the 
absence of a mass term, i.e. in the ensemble of fixed applied 
force /. It is based on the assumption that, as the force / 
is increased towards the threshold fc, the mean number of 
avalanches per interval of force increment is not singular. 
Translated into the present setting it means that 



rn^dw 



Po 



(160) 



remains finite as to ^ 0. (We use the notations of Section 
|III C| where N^^ is the number of shocks and po the shock 
number density). Note however that in the present setting, 
although the increase of the force is iv?dw, it does not bring 
the system closer to criticality, which is achieved by decreas- 
ing the mass. This is why we had at the start two independent 
exponents p and r. The STS symmetry provides a first re- 



lation (23 I while the conjecture of the finiteness of ( 160 1 is 



equivalent, using p9] l, to the conjecture p = 2 (at least for 
T > 1). Using (29 1 again, this implies the value r = r^. 



Finally note that a similar conjecture was put forward in the 
context of sandpiles 1 18 19|, and if one admits the connection 
to the depinning of a periodic interface Ii24ll75l it is equivalent 
to the conjecture t — t^^q for the RP class. 

As discussed in Ref. Il67l . we expect the result for r for 
statics (shocks) and for depinning (avalanches) to coincide to 
one loop. Hence our result indicates that both indeed coincide 
with the conjecture r = r.^ to one loop. What happens be- 
yond one loop accuracy remains to be elucidated. We know 
that statics and depinning differ at two loop, e.g. for the C, ex- 
ponent. There are various possibilities for either the statics or 
the depinning to agree, or disagree with the conjecture at two 
loops and for c? = 2, 3: Since (^cp > Cstat, one possibility is 
akeady excluded, that both conjectures are true and that the 
exponent r is the same in both problems. Some claims that 
static and dynamic avalanches belong to the same universal- 
ity classes were indeed put forward in |22| for the RF Ising 
model. Finally note that our result can also be written as 



and a remaining challenge is to compute the higher orders. 
Until now we have failed to find, within the field theory, a 
symmetry or a mechanism which would imply p = 2 to all 
orders. 

Note that exact results for d = depinning ll86ll . i.e. of a 
particle, yield t — for the so-called Gumbel class which has 
C — 2, hence agreement with the conjecture up to logarithmic 
corrections present in that case. However, the result for the 
other classes (Weibul and Frechet) have also r = but C 7^ 2, 
hence violating the conjecture in its present form. For the 
statics, exact result for the d = limit of the RF class, the so- 
called toy or Sinai Brownian energy landscape, has t = 1/2 
and — A/3 which again satisfy the conjecture |,62J . 

To summarize, for the cases of interest, our result is, with 
e — 4 — d and up to 0{e^) corrections: 



1 

12' 



- < 0.0989627 e 



1 



RF 
RB 
CDW 



(162) 



where we have used the 1-loop values for ( in the statics, re- 
called in Appendix [B] To 0(e) these are the same as for de- 
pinning, except that in that case the RB class does not exist, 
since 2-loop corrections make it flow to RF. 

It is useful to quote the expected result if the conjecture 
holds. In d = 3: 



r = < 



In d = 2: 



T = < 



7 

5 

1.37783 ±0.00005 
4 
3 



5 

4 

1.17987 ±0.00037 



RF 
RB 
CDW 

RF 
RB 
CDW 



(163) 



(164) 



In d = 1: 



RF 
RB 

CDW 



(165) 



p = 2 + 0(e2) 



(161) 



For the RB class we have used the 2-loop result ( — 
0.208298042e + 0.00685826^ + 0{e^) and constructed three 
different Pade appro ximants, using the constraint that Cd=i — 
2/3. This gives (^=3 = 0.21454 ± 0.00028 and Cd=2 = 
0.43864 ± 0.0011. The error bars denote the la spread of 
the three Pade approximants. For RF we use — (4 — d)/3 
and for RP C = 0, hence there is no uncertainty. These are the 
values for the statics. For depinning, if the conjecture holds, 
the values for the RP (CDW) class are the same as given here. 
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For non-periodic interfaces C = 56(1 + 0.14331e) + O(e^) 
and one gets r^^s = 1-409 ± 0.001, rd=2 = 1-31 ± 0.01 
and Trf^i — 1.2 ± 0.1 (in the latter inserting the numerically 
measured value for ( = 1.25 yields instead Td=i — 1.11). 



order by order, by first inserting in the term proportional to a 
the MF solution, i.e. the solution of ( |150[ ) at a = 0. Solving 
for Z then yields 



B. Distribution of avalanclie sizes 

We now compute the avalanche-size distribution. To do so, 
we have to invert ( |150| l. To first order in a this can be done 



^ 1 . _ . ((3A + yr^-l)log(1^2(2A + yr^-l)). ^ 



(166) 



Expanding in A gives 

Z(A) = (A + A^ + 2A3 + 5A4 + 14A5 + 42A6 + O (A^))+ (^2A3 + ^ + 54A5 + "^I^ + q (a^)^ ^+0 (a^) , (167) 



consistent with our previous result ( 144 1 if one uses the values of the i„ given above. The challenge is to find a distribution P(JS) 
;nerati 

Z(A) 



which generates all terms in ( 166 1. By trial and error, one arrives at the following integral representation: 
1 



AS (e^^ - 1) 5-3/2 

1 + ^ (log(S')S' + 75 + 45' - 87^75 - 6 log(5) - 67 + 4) a + O (a^ 



(168) 



To leading order in a, using the definition ( 152 1, this yields 
our final result for the avalanche-size distribution at one loop: 

|5\ 



SSI 
2^ 



S:^-'^AS~^ exp C 



B 

J 



S 



(169) 



for 5 > 5,„in 

A = 1 + ^(2-37E)a 



with coefficients: 
1 



B = l-a 1 



7e 
4 



C = 



-3(1-C.)e 



(170) 



7e — 0.577216, and the exponents read: 



3 
2 

1 - 



1 



(l-Ci)e- 
(l-Ci)e. 



(171) 
(172) 



a ^ 1 
4-1+12^ 

Note that the decay of large avalanches becomes stretched ex- 
ponential, faster than exponential, with an exponent S > 1. 
For the RF class Ci = 1/3, hence ^ = 1 + e/18 + 0{e^). 
This goes in the right direction to match the exact result S — 3 
for the RF class in rf = 0, but the magnitude seems under- 
estimated by the 1-loop formula, which suggests large higher- 
loop corrections for this exponent. Note the pronounced bump 
in the plot of S'^P{S) on figure |2] This bump is indeed seen 
in numerical simulations 1261 . 



C. Normalization and scales in tlie distribution of sizes 



One can check that the final formula ( 169 1 is properly nor- 
malized, i.e. dSSP{S) = (5), to leading and first order 
in e. Note that since t > 1, this formula does not give in- 
formation about typical avalanches, which are of the order of 
the cutoff 5min, but about larger ones which control the mo- 
ments (S^) with p > pc — T — 1. It is then universal, i.e. it 
does not depend on the details of the small-scale cutoff Smin- 
However the average avalanche size (5), which appears as a 
factor in the distribution ( |169| l, is non-universal and cannot be 
computed from this theory. In Section V F we showed that it 
behaves as (S) ~ S^^S^^ but this is all one knows about 
it. In comparing with numerics or experiments, one replaces 
(5) by its actual measured value. Then there is only one free 
parameter left, the global scale Sm- Note that this scale is 
given by Sm — {S'^) /2(5), hence it can also easily be mea- 
sured. We recall that its value is predicted by field theory via 
the exact relation 



Sm — 



|A'(0+)| _ K'^/^ |A'(0+)| 



(173) 



el2 = 2(47r)^''/2p(3 _ ^) is a d-dependent number A ques- 
tion is how predictive this formula, and how universal this 
scale is. The answer depends on the universality class, RP, 
RForRB. 
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For the RP class, C = and the fixed point for |A'(0+)|/a, 
where a is the period, is a universal number. To two loops, 
this number is equal in the statics and at depinning and reads 
e/6 + e^/9 + 0(e'^). In general K and m will experience small 
corrections for any e.g. a lattice model which does not exactly 
satisfy the STS symmetry. They can however be measured 
from large scale measurements on the system. Once they are 
extracted, then Sm can be predicted. Alternatively, one can 
construct ratios which, for the RP class, are universal, for in- 
stance. 



aSjn 



a|A'(0+)| _ a|A'(0+)| 



A(0) 



A(0) 



= 6 , (174) 



where u is the center of mass and v? the variance of its fluc- 
tuations. Its definition is u'^L'^ = u{Q)u{x). Note that in 
the statics there is no higher order correction, since the fixed 
point for R{u) — a + h[u{l — u)Y with some constants a and 
h which drop out from ( |174[ l. Replacing by the depinning val- 
ues gives r = 6 + 2e + 0(e ). This universal ratio can thus be 
used to distinguish statics from depinning. 

For the RF class, = e/3 in the statics, as recalled in Ap- 
pendix |B] The fixed point values contain a scale ^ (in the di- 
rection of u), which can be fixed if one knows the amplitude 
of the random field cr = /q dw A (u) . ct can be retrieved from 
large distances in experiments, and is a parameter of simula- 
tions. There are various amplitude combinations which can 
be studied, depending on whether one is willing to measure 
TO and K, see Appendix[B]and Ref. ||87J . The nicest universal 
ratio is: 



TKF 



a|A'(0+)| 



= 71 



2 
9 



A(0)2 
71 +72 ) +0(e^) 



(175) 



with 71 = 0.775304245188 and 72 = -0.13945524, where 
one can check from Appendix [B] that all dependence on to, 
K and the scale ^ cancel (thanks in part to the strict equality 
C = e/3). 



Finally, for the RB class (and for the RF class for depin- 
ning) there is the least universality: | A'(0+) | is non-universal, 
hence measuring K and to will not be enough. One has, for 
instance, also to measure vP'L'^ to obtain the combination: 



A'(0+)2 m-'^K'^''^ [A'(0+)]^ 



u-^Ld to4A(0) 
with the universal ratio 



e/2 



A(0) 



(176) 



A'(0+)^ 
A(0) 



= 0. 5834056 + 0.2942056^ + 0{e^) . (177) 



Of course, in all these cases there are other possibilities for 
interesting universal ratios. If one measures for instance the 
function A(u) as in Ref. [26 1, one can in some cases get rid 
of measuring K and to. Eventually, a direct numerical test 
of the relation between Sm and the cusp would also be wel- 
come. Analytical solution of toy models in d = L86j has 
also successfully tested this relation. 



D. Moments and universal ratios 

Having obtained P{S) we can compute its moments. This 
is useful for comparison with numerics. Direct integration of 
equation ( |171| i gives 

-V{n-\) - 5) + log(4) + 7) (n - 2) + 4n] } 

+ O (e^) (178) 

a formula valid for any fixed real n > 1/2. From this we can 
extract the universal dimensionless ratios: 



(179) 



(5") 2 2n-3 
e nr(n-|) + V^r(n-l) 

3^^-^^) (2„-3)2r(n-|) 
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for any real n > 3/2, with Ci = 1/3 for RF, Ci = for RP 
and Ci = 0.208298042 for RB. The lowest-order integer ones 
are: 



IX. SPATIAL STRUCTURE OF AVALANCHE 
DISTRIBUTIONS AND SELF-CONSISTENT EQUATION AT 
NON-ZERO MOMENTUM 



r2 = 3 - e(l - Ci) 

5 5 , 
"^^3-27^(^-^^)- 



(180) 
(181) 



Another useful form when comparing to numerics for n 
near r, is to isolate the simple pole divergence which occurs 
in any dimension: 



r,,, = h B„,d 



n — T 



(182) 



12 



(183) 



Br,M = 1 + e(l - Ci) 



7rr[n- 1/2] - V^r(n- 1) 



Note that for comparison with numerics it is useful to es- 
timate the corrections due to the small scale cutoff, assuming 
it just cuts the previous result at S„iin- The correction to the 

^ are thus 



{S}Sl 



above result for the dimensionless quantity 

of order {Smin / S„i)"^'^^^ / {n — t + 1) and an approximation 
to r„ is thus: 



A, 



(184) 



For the RF class one finds either from a direct expan- 
sion, or using the (0, 1)-Pade: = {2.333,2.45455}, 
ra = {1.5432,1.55172} (d = 3), r2 = {1.667,2.07692}, 
rg = {1.41975,1.45161} (d = 2); ra = {1,1.8} and 
ra = {1.2963, 1.36364} (d = 1) and ra = {0.3333, 1.58824} 
and ra = {1.17284, 1.28571} {d = 0). 

This can be compared with the exact result ll62l for RF dis- 
order ind — 0: 



r2 = 1.2978 
rg = 1.17776 



(185) 
(186) 



Using the exact result as a constraint yields three Fade aprox- 
imants and the results: 



and 



r2 = 2.404 ±0.009 d=3 
r2 = 1.935 ±0.021 d^2 
r2 = 1.571 ±0.022 d=l 



ra = 1.5427 ±0.0012 d=3 
ra = 1.419 ±0.003 d = 2 
ra = 1.297 ±0.003 d=l 



(187) 
(188) 
(189) 



(190) 
(191) 
(192) 



where averages are over the three Fade approximants and error 
bars are the corresponding one sigma deviations. 



In this section we introduce generating functions which en- 
code for the spatial correlations in the avalanches. An explicit 
calculation is performed at the level of the improved tree ap- 
proximation (mean field). It exhibits an interesting connection 
to instanton calculations in a cubic field theory. 



A. Generating function 

To obtain information about the structure of avalanches in 
internal space one may define for avalanche i 



Sf 



(193) 



where (j){x) is a given function. One recovers the standard 
definition of size for 0(a;) — 1, i.e. Si — Sj. One would like 
to compute averages such as 



1— y (e^/. ~ l)S{w - w,) 

-7)1 T-l ^ — ^ 



1 

L 't'(^) Jo 



.AS* 



l)p^{S^), (194) 



where Pcf,{s) — poP^{s) is the density of avalanches with 
= s, and Pcf,{s) the normalized distribution of 5*0. Note 
that STS implies (l){x){uj;{w) — w) = for any 0(a;) and 
w. Hence, taking a derivative w.rt w one obtains 



(195) 

Note that po, the total density of avalanches is independent of 
(p, hence one also has the exact relation 

(5^) = (5) L-^ I 0(x) . (196) 

Extending the arguments leading to ( |A4[ ) and ( |A5| l we write 

dS'"^P0(S"^)(e^'^*-l-AS'"^). (197) 
Hence one needs to compute the generating function 

G'^(X) = - qHI^ 'l'ix)i^:^iw)-w-u^{o))) (198) 

(199) 
(200) 



G"^(A) = Z'''{\)w + 0{w'^) 
^nA) = ^((e^^*)-l-A(5^)) 



for w > 0, with again Z"^(A) = A + Z'^(A). Note that here 
we consider only a uniform — w. 
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B. The self-consistent equation 



Multiplying with Y' and integrating once gives 



Let us now study this quantity in the improved tree approx- 
imation (also called mean-field above). We can show that it is 
given by 



1 



where Z!^{X) is solution to the self-consistent equation 



(201) 



Z^{X) = A</.(x) + |A'(0+)| / g,.yg,..Z^{\)Zt{\) 



(202) 

<7x = Jf. gk&^^^ is the free (elastic) propagator. One can check 
that for (t){x) = 1, Zx{X) = Z{X) and one recovers the tree- 
level recursion given in the text for Z{X), and, upon replacing 
A'(0+) 1 and gk — > gk — gk\m=i, the one for the rescaled 
function Z{X). A derivation, resumming all diagrams, is given 
in appendix |F 1| 



We now go to rescaled quantities, using (91 1, i.e. Z(X) 



Z{XSm)/ S„i- Then we can write for an arbitrary function 
(j){k) in Fourier space: 



2~~i 2T1 2ZqWZk-qW 

q + (fc — q)^ + 

(203) 



Zfc(A) = xm 

In real space, this is 



Z,iX) = X(b{x) + / g,^^g,^yZ,iX)Zy{X) . (204) 

Jxy 

We note that this self-consistent equation is simplified by 
defining Yfc(A) :— m^gkZk{X), or in real space 

Y{zm, A) f 5.-.^x(A) (205) 

J X 

ip{mz) 0(z) . (206) 

This results in (suppressing from now on the explicit depen- 
dence on A when convenient) 



(-V^ + 1) Y{x) = Xif{x) + Y{xf 



(207) 



C. Solution for </? localized on a codimension one hyper plane 

The function Z{X) will be qualitatively different, depend- 
ing on whether is extended on the scale of the inverse mass 
(absorbed in x), or is a (5-distribution. We now study one spe- 
cial case, namely f{x) — 8{x\). The function Y{x) will then 
be constant along the directions a;2, . . . , x^, and for simplicity 
of notation we will denote x\ x, and suppress Xi, i > 1. 
Thus we effectively consider a 1 -dimensional problem. 



Eq. (207 1 can then be integrated analytically. Consider first 



the homogenous equation (A = 0) 

Y"{x) = Y{x) - 



(208) 



^\x)f = const - -Y{xf + Y{xf 
3 



(209) 



If ( |208| l is viewed as the equation of motion of a particle, 
then setting const will be the solution which has zero 
kinetic energy at the saddle-point F = of the potential 
= IF'^ — iy^. It is the unique solution which de- 
cays (exponentially fast) to F = for x ±00. The other 
solutions are either oscillating or unbounded. Integrating once 

more for const = gives —2 arctanh 

alently 



\/3 



x, or equiv- 



Y{x) = Y^{x) ■-- 



cosh(x) + 1 



(210) 



where the center of the solution has been chosen to be at a; = 
0. 

A symmetric solution Y(x) — Y(—x) of (207 1 with 
^{x) — 5{x) can now be constructed as follows: 



r(x,A) = 



for a; > 



cosh(a:: + a;o(A)) + 1 
Inserting into ( |207| i and integrating from —8\a 8 gives 

- 2d&Y{8, A) = A + 0((5) . 
Thus, in the limit of (5 — > 0, 

sinh(a;o(A)) A 



(l + cosh(xo(A)))2 



6 



(211) 



(212) 



(213) 



Let us consider the unique real branch, s.t. xo(A) 00 when 
A ^ 0, which gives Z'^(A) ^ in the same limit. Increasing 
A from 0, the solution breaks down, when A reaches Ac with 



Ac 



2 



(214) 



such that cosh(a;o(Ac)) — 2. Now we need Z"*(A) defined m 
, or the rescaled version Z'^(A) 



(201 



Z^(A) = 



= j^^jyV^ + l)Y{x,X), (215) 
where from the first to the second line we switched to dimen- 



sionless variables. Inserting (207 1 and using that f{x) — 6{x) 
yields 



Z'^iX) = f X^{x)+Y^{x,X) 

J X 



= A + 
= A + 12 



da; 
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[l + cosh(a; + a:o(A))]2 
1 + 3e^o(A) 



[l + e=^o(A)]s 



12 



1 -(_ qXo(X) 



(216) 
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I 



10 



FIG. 3: Z'^{X) for 4>(x) — 5{x), as explained in the main text. 



where in the last line (213 i was used. Therefor e xa( X) can be 
expressed in terms of Z'^(A), and inserted into (213 1, with the 
result 



A = 



Z'I'iZ'l' -6){Z'f' -12) 



72 



(217) 



Note that this result was derived for A e [0, A^], Z"^ G [0, 6 - 
2\/3], which corresponds to the red (solid) branch in figure 
[5] and then (implicitly) continued analytically to all A and all 
Z^. We note that the green part of the curve (long dashes) 
can also be obtained analytically, using the second point for 



which (213 I is satisfied, and the two solutions for negative A. 



One easily checks that the symmetry Y{x) — Y{~x) implies 
point-reflection symmetry around the point (A, Z'^) ~ (0, 6). 
This already strongly suggests that A is a third order poly- 
nomial in Z'^, symmetric around the point Z'^ = 6 (in our 
chosen units). 

For large negative A, |Z*(A)| - y|A|, thus fr om our dis- 
cussion following ( 153 1, Pc ~ 5, and the avalanche exponent 
T is 

r*=l- (218) 

Note that this value is consistent with a generalized conjecture 
that we can put forward, r"^ = 2 — where = d ~ 1 
for a codimension one subspace, inserting d = 4 and C = 
since the above result was derived at mean-field level. Again, 
much work remains to validate or invalidate this generalized 
conjecture (e.g. a two loop calculation). 

Interestingly, the probability distribution is non-trivial and 
different from the standard mean field one (i.e for (p — 1). It 
takes the form, for 5"^ > St' 



1 



i-) 



(219) 



Here the avalanche sizes are defined as Sf 

— f SI. n- ^ ~ m""^"^ and Sm is the same 



quantity as defined in (93 1. The distribution p'^(s) can be 
calculated as follows: 



y(s) 



dA 
27ri' 



-As 



Z*{\) 



dZ dX{Z).,z 



(Z)s, 
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2m dZ 



dx {3x^ 



X [a:sin(3sa;(a;^ + 1)) - cos(3s2;(a;^ + 1))] 



2Ki 



(220) 

77 S 

Here s = S't'/S^ and (s) = (S"^) /S*™, and Ki is a Bessel 
function. The steps of the derivation are: a change of variables 
from A to Z = Z'^(A); a change of variables Z = 6i{x — 
i) and combining the integrand for x and —x. For large s 
the asymptotic behavior is similar to the standard mean-field 
result P{s) ^ s' 



-3/2, 



'l\5/2 



(221) 

This was expected, since in both cases the solution for Z{\) 
ends at a A^ with a square-root singularity. For small s the 
asymptotic behavior is in accordance with ( |218| i: 



{s)-'p'^{s) 



,2/3 



7rs4/3 ^7rs2/3 27r 
Finally, from \22Q\ one obtains the moments: 



3«/2p I 



6) 

2n 



(222) 



(223) 



for n > 1/3, from which the universal ratios r„ can be com- 



puted. In particular that ((S"^)^) / (5^) 
2S'„ for 6^1. 



hSm instead of 



Note that for cj) localized on a hyperplane, ( 196 1 implies 



Lm 



(S) 



(224) 



This is understood from the observation, that only a fraction 
1 /{Lm) of all avalanches leads to an advance of the interface 
(avalanche) in the hyperplane. Thus, if one were to define 
an avalanche distribution by considering only the hyperplane, 
one would naturally choose a different normalization. 



D. Diagrammatic expansion 

Here we have solved the codimension one case, i.e. d^ = 
d—1. More generally one may consider the case d^ = d — d' , 
i.e. a function (^(x) localized on a rf' < d dimensional hyper- 
plane. This amounts to study Eq. (207 1 in dimension d' . This 
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can be done graphically, as (207 1 can be used to define, in any 
dimension d! , a diagrammatic expansion: 



Z{\) = A + A 




(225) 



where independent d'-dimensional momenta flow in each 
loop, the solid lines are propagators — + 1) and 

each loop contains a momentum integration. Note that this 
expansion is valid, i.e. IR- and UV-convergent, in all dimen- 
sions d! < A (hence any d > for the original problem). In 
d' = it reproduces (89i. In d' = 1 the integrals are eas- 
ily computed by switching from Fourier to real space where 
gx = ^e^l^l, and one obtains 



(226) 



This is in agreement with the analytical result (217 1 

A2 A3 SSA"* A^ lOOlA^ 4A^ 
T ^ 576 27 ^ 41472 ^ 243 

46189A8 55A9 5311735Aio 

(227) 



Z(A) = A+- + - + 



3981312 6561 



859963392 



It is amazing that the quite non-trivi al ser ies ( |225[ ) is indeed 
resummed by the instanton solution, (21 1 1 in the case d! = 1. 



Study of higher d' and extension to the loop expansion is left 
for the future. 



X. NON-LOCAL ELASTICITY AND AVALANCHE 
DISTRIBUTIONS FOR A CONTACT LINE 

The calculations performed in this paper can be extended to 
a broader class of elasticity functions in One possible 
generalization is 



(228) 



where 7 = 2 corresponds to the choice studied above of local 
elasticity, while the case 7 = 1 is relevant for long-range elas- 
ticity as in contact lines of fluids. In fact a more realistic form 
in that case is 



(229) 



The mass is provided by capillarity, i.e. the interplay of 
surface tension and gravity. We keep the notation = 
.95=0 = A*^ for the curvature of the quadratic well, while dis- 
tinguishing it from fj, which is an inverse characteristic inter- 
nal length in the interface (the inverse capillary length). One 
parametrization which contains these three cases is 



9q ^ = m'^5(9/m) 



(230) 



for some function g{y) ^ y'^ at large y, and g{Q) = 1. Of 
course ot her sc ales may be present in a realistic problem, and 
the form (230i assumes that fjT^ is the single largest internal 



length scale which cuts off the avalanches. 

The upper critical dimension dc in all cases is given by the 
UV divergence of 



(231) 



where e — dc ~ d, dc — 27. The constant Cd^j = e/2 with 
h = 5gU=i remains finite at 4. 

With these choices all calculations of this paper are easily 
extended. Appendix [B] can still be used, up to trivial changes, 
e.g. replacing to by /i in all final formula there, including in 
the definition ( |B3[ ) of the rescaled disorder (we can set K — 
1). One has again 



-A"(0^ 



-(l-Ci) + 0(e') 



(232) 



to one loop, with the same values for (i, independent of 7 
(and g{y)), and of course now e = dc — d everywhere. For 
more details, including the calculation of the exponent ( (and 
z for the depinning) to two loops for the non-local elasticity 
see Ref. 13611491 . Here we compute the avalanche distribution 
to one loop. 



To do that one starts again from Eq. ( 151 1. It generalizes as 
follows: 



Z(A) =: Z= X 



Z2 



(233) 



Z 



eh Jk {~gl' - 2ZY ~gl' - 2Z 



^-Z~gk- 'iZ'gi 



(234) 



where g^ = gk\^l=l■ Note that the definition of the rescaled Z 
is still given by \\52\ with now 

TO-*A'(0+) = 



5, 



-27A'(o^ 



(235) 



The calculations are performed in AppendixjE] For the choice 
of elasticity ( |228[) we find again, for any 7, the distribution 
P{S) given by (68 1, with amplitudes 



A = 1 + -^(2-37E)a 
47 

C — \/TTa , 

7 



B = 1 



27 



and exponents 

3 

"^2 



1 

47 



(236) 



(l-Ci)e + 0(e2) (237) 



(5=1 



1 



l+^(l-Ci)e + 0(e^) 



(238) 



47 
a 

27 ' 67 ' 

Of course everywhere e — 2j — d. Note that the 1-loop result 
for the avalanche exponent t is now compatible to 0{e) with 
the generalized conjecture 



7 



d + C 



(239) 
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This conjecture can again be reexpressed as 



given that the scaling form ( 22 1 generahzes into 



p{S) = i'^m" 5-^(5'/+'^) , 



(240) 



(241) 



and that the STS relation (21 1 still holds, hence now the expo- 
nent relation (always true) reads: p = 2(2 — r)((i + C)/l- 

In the case 7=1 and for the form ^229\ more suitable to 
describe the contact line, one finds again the same exponent t 



(237 1, however the shape of the distribution is different. We 
find, see Appendix IE) 



P{S) 




X exp C" 



(242) 



for S S'min, with amphtudes 



A' = l + i(8-37E)a 



B' - 1 - 2«(7E - 2) 



C 



3 r 



The exponent 5 is now given by 

<5' = l-^ = l+i(l-Ci)e + 0(6^) 



(243) 



(244) 



Note that the presence of the constant D' suggests that the 
pre-exponential power law at large S is not S^"^ any more. 



XI. CONCLUSION 

In conclusion we have succeeded in computing from the 
FRG the distribution of shock or jump sizes which occur in the 
minimal energy configuration of an interface pinned in a ran- 
dom potential and tied to a spring of varying position. These 
are the static analog of the avalanches which occur in the dy- 
namics if the interface is instead pulled from a metastable con- 
figuration to the next one. Hence it opens the way to the same 
calculation in the dynamics, performed in ||67l , which yields 
very similar results. Shocks in the statics are often called static 
avalanches as there are many analogies, as well as some dif- 
ferences. 

We believe that this is an important achievement. First be- 
cause the FRG has been around for a while but it was not 
understood previously how to extract the avalanche statistics 
in a controlled way. In the sandpile literature this is still an 
open question despite many exact results for other quanti- 
ties. It turns out to be conceptually simple (a posteriori) to 
extract these distributions from the FRG. In fact the beauti- 



fully simple relation ( 36 1 between the cusp of the FRG func- 
tion A(it) and the avalanche-size's second moment, unveiled 
in this work, gives a very transparent physical picture of the 



cusp. Similar relations hold for all moments, and the chal- 
lenge is to sum them. This part is a priori technically difficult 
but some surprising simplifications occur in the calculation, 
which in the end lead to a simple self-consistent equation, 
with an interpretation in graph theory. This suggests a, yet to 
be discovered, simpler and presumably more powerful struc- 
ture behind the present state of the art of the theory. 

The 1-loop result for the probability distribution is equally 
striking. First the avalanche exponent r is found equal to or- 
der 0(e) to a conjectured form, put forward for the depinning 
transition. Note that our result is a first-principle derivation 
of this exponent. Hence the conjecture is confirmed to one 
loop. It is then of high interest to look for possible deviations 
to two loop, both in the statics and dynamics. Alternatively, 
if the conjecture is true in general, it would be interesting to 
derive this from first principles in the field theory. Besides the 
exponent, we obtained the general formula for the distribution 
P{S). This includes the precise way in which it is cut off for 
large avalanches, at a cutoff scale Sm- These large and rare 
events are the important ones in terms of moments of distribu- 
tions, and, in real life, if one is interested in e.g. earthquakes. 
We predict the scale Sm and the distribution, which is insen- 
sitive to details of the model at short scales (up to a single and 
measurable global factor). Of course the distribution depends 
on the details of the model at large scales, here we use mainly 
a quadratic well to cut the size off, but this dependence can 
explicitly be computed within the theory. We illustrate this 
point by computing it for contact-line depinning which has a 
more complicated elastic energy. 

Some progress was made to study the spatial structure of 
avalanches. Additional definitions of local avalanche sizes 
was given which integrate information about jumps within a 
subspace of the interface. It was shown that to obtain their 
probability distribution at mean field level one must solve an 
instanton problem in a cubic field theory. An explicit solution 
was found for a subspace of co-dimension one, leading to a 
novel exponent r = 4/3 at mean-field level, i.e. in d = 4, and 
a novel size distribution involving a Bessel function. It is dif- 
ferent from the usual mean field size distributions, recovered 
here, with r 3/2, and an exponential function. 

There are many open interesting questions which can now 
be addressed, both within the statics and the dynamics. First, 
one would like to know more about the spatial structure of the 
avalanches, and their correlations. A first step is to extend the 
instanton calculation to regions of arbitrary shapes and co- 
dimension, and to study it within the loop expansion. The 
problem can also be extended to a manifold with N compo- 
nents. Consequences for hysteresis loops of magnets deserve 
a new study. Applications to earthquakes are also of interest, 
especially if avalanche correlations can be handled (pre- and 
after-shocks). Plastic avalanches can now be studied, e.g. in 
the framework of Ref. 
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APPENDIX A: SHOCK-SIZE DISTRIBUTIONS 

1. Generating function 

We start again from the shock decomposition Eq. ( [T6| . The 
field e^^ lu(w)~w] jj^g jump discontinuities at each shock po- 



sition, where its value is multiplied by a factor e 
one can write 

i 



Hence 

(Al) 
(A2) 



where LMv^H = E,<.Sj and L^O = E^<^S,. 
The labels of the shocks are ordered according to thek spa- 
tial position Wi. Let us multiply these equations by the field 
gAL [«(tu)-«;] i-^]jgjj ^j- ^ different point, i.e. consider 



-XL''[u{w2)-W2] 



(A3) 



where wi and W2 are generic points and different from the 
Wi (even though the notation might suggest otherwise). If we 
consider the limit W2 = in < A3 1 and average over disorder 



the term containing u{w^ ) cancels and one obtains: 



lim {d^,^ + XL"^) e>-L''[u{w,)-w,-(u{w2)-w2)] 



i 

= I dSp{S){e'' - 1) . 



(A4) 



A similar equation can be derived starting from ( A2 1 and con 



sidering the limit W2 — ■ Using translational invariance 
one finally obtains 



w=0+ 



a„e^-^''["('")-"'-"(o))] 



d5'p(S')(e^^ -1- AS*) 

(A5) 

J dSp{S){l - e-^^ - XS) , 
(A6) 



where we have used (20i. One thus recovers Eq. (44i in the 
text. Another method to derive the part proportional to \w\ is 



to multiply ( Al i at two different points (and setting A ^ — A 



for the second) and extract the single shock contribution: 



(9^2 + XL'^){d^^ + AL"*) c>^L''[u{wi)-w^-(u(w2)-w2)] 



= 5{wi — W2)^{e^^^ - - l)5{wi — Wi) 

i 

+ smooth . (A7) 
This yields, using translational invariance: 



{~dl + X^ L^"^) e>^L'i{n{w)-tu-u{o))) 

= 25{w) J dSp{S)[l~ cosh{XS)] + smooth . (A8) 



Integrating twice one recovers Eq. (44 1, but only the coeffi 



cient of |w| is determined. To determine the coefficient of w 
one needs another equation, e.g. as above, or has to use the 
symmetry under w ^ —w. 

Note that this can be generalized as follows: Consider a set 
of Xk with = 0, then one has: 



fc=i 



S{'wi — W2) ■ ■ ■ S{'wi — w„) 

/n 
dSp{S) Yl {e^"^ - 1) + less singular (A9) 



fc=i 



2. Multi-sliock size distribution 



The multi-shock size distributions are encoded in the higher 



pinning force cumulants ( 13 1. For instance the third cumulant 
can be written as: 



: m^L^'^{u'{wi) - l){u'{w2) - l)(u'(w3) - 1)' 
m^L-'' ( dSS'^p{S)5{wi - W2)5{wi - w^) 



(AlO) 



j dSidS3SiS3Pc{Si,S3,wi,W3)S{wi - W2) 



2 perm 



+ m'^L J dSidS2dS3SiS2S3Pc{Si, S2, S3, wi,W2, W3) 

Pc{Si,...,Sn,Wi,...,Wn) = ^ 



all distinct 



(All) 

which is easily generalized to any cumulant. Note that the 



connected parts and the — 1 substractions conspire, using ( 20 1, 
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to give the correct connected shock distributions (such that 
scaHng with volume is straighforward). Such formula show 
explicitly the structure of the sigularities expected in the force 
ri-cumulants as p points are brought together, as a conse- 
quence of the assumption of a finite density of shocks (dilute 



shocks). In principle it allows to check this assumption by 
computing all moments within the e expansion. 

The generating function generalizing Z, and allowing to ex- 
tract pc{Si, S2,wi,'W2), which is a function of wi — W2, is 



constructed in analogy with (44 1 and (45 i: 



G{\i,\2,Wi,W2,5i,52) = Z{\i,\2,Wi - W2)5i52 + ■ ■■ (A12) 
G{\i,\2,Wi,W2,5i,52) = L-'^(e^iL''[u(w^+S^)-u{w^)-S^] „ 1) (^q\2L'^[u{w2+5^)-u{w^)--5^\ _ l^^" (^O) 



for(5i,(52 > 



r 



APPENDIX B: REVIEW OF BASIC FRG RESULTS 

We review basic equations for functional RG, mostly at 1 - 
loop order, making explicit the dependence on the elastic co- 
efficient K, which usually is set to one. Details (mostly at 
if = 1) can be found in ||36l|49l. The FRG equation for the 
function A(m) to one loop is: 



mdmA{u) 



1 



2-\rr 



(Bl) 



where /„ = + m^) with an impUcit UV cutoff A. 

One has —mdmh = ^m^I^- One should distinguish d < 6 
and d > 6. The latter is dominated by the UV cutoff and is 
briefly discussed at the end. For d < 6, in the limit m <C 
A one has h = hm<^-^ K~<^/'^ with /„ = /^^(fc^ + l)-" 
(defined with infinite UV cutoff whenever convergent). Thus 
h = (ei^2)/4 with 12 = (47r)-''/2r(2 - f ). Hence for d < 6: 



- mdr,j2 = 4.m^h = m'^-'^{ei2)K-'^'^ 



(B2) 



where the combination e/2 = 2(47r)^''/^r(3 - |) is wefl de- 
fined for all d < 6 (the pole at d = 4 is suppressed). One 
defines the rescaled (dimensionless) function A(u) through 



A{u) = — ^m""^^A(Mm^) 
eh 



(B3) 



with e/2 = l/(87r^) in d = 4. It satisfies the dimensionless 
FRG equation: 



md^A{u) = (e - 2C)A(u) + CuA'{u) 
-i[(AH-A(0))2l 



(B4) 



This equation is valid for all d < 6 and yields a fixed point 
A*{u) = 0(e) for d < 4 with ( — Cie and a few universality 
classes. The two loop equation, not reproduced here was also 
analyzed. We recall the main results: 

(a) random-bond class: 



C = 0.208298042e 
[-A*^(0+)]^ 
A*(0) 



O.OOeSSSe^ + 0{e^) 
0.583405e + 0.2942056^ + 0{e^) 



(B5) 
(B6) 



Since (B4i is invariant under A{u) ^^^A(^u), A{u) con- 
tains one non-universal scale, and only the above ratio is uni- 
versal. 

(b) random-field class, with R{u) — —cr\u\ at large u: 



c 



(B7) 



A*{o) = -e 



-A*'(0+) = -e 



where ^ is the non-universal scale introduced above. Integrat- 
ing A{u) from to cx) yields gU, Eq. (4.39): 

K-''^^{ei2)a = [71 + £72 + 0{e^)] (B8) 

with 71 = /g dy v^2(y- 1 - Iny) = 0.775304245188, 
72 = -0.13945524 (see |49|, Eq. (4.41)). 

(c) periodic class with period u = a: 



C = o 

A*(0) ^ 
-A*'(0+) = 



e e 
6+9 



For all classes: 



A*"(0+) = e(l-Ci)/3 + 0(e2) 
In dimension d = 4 one defines: 

A{u) = A(uln(mo/m)"^i)(ln(mo/TO))"i+^'^i 
and setting C = and e = in (|B4|i one finds: 



(B9) 
(BIO) 

(BID 

(B12) 
(B13) 



edeA{u) = (1 - 2Ci)A(u) + CimA'(u) 
-M(AH-A(0))2T 



(B14) 



where £ = ln{mo/m). Hence A{u) converges to the same 
fixed point A*{u) — A*{u)/e as for e > 0. This yields for 
the original function A{u) in d 4: 

A{u) = 87r2i^2A*(uln(mo/TO)~''^0(ln(TOo/™))"^+^''' 

(B15) 
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up to subdominant terms of the order of 1/ ln(TOo /"^). Here 
1/mo is a non-universal scale, presumably (at least) of the 
order of the Larkin scale. In dimension 4 < d < 6, we set 
again C = 0- The equation for A"(0) reads: 

-m9„,A"(0) = -{d - 4)A"(0) - 3A"(0)2 (B16) 

for an analytic disorder (A"(0) < 0). Hence if the bare dis- 
order is smooth one needs |Aq(0)| > (d — 4)/3 to generate 
a cusp and metastability. Hence for sufficiently strong bare 
smooth disorder, or with (even weak) rough bare disorder, a 
cusp is generated. Eventually as m ^ the flows converges 
back to the attractive fixed point A(u) — as: 



AH - ( £^ ) K*{u) 



A{u) 



m 

mo J 



(B17) 
(B18) 

(B19) 



where now A* (u) is non-universal and depends on details of 
the FRG flow at intermediate stages. It has a cusp, unless one 
starts from smooth weak disorder. Note that for d > 4, if 
we make the natural assumption (as for d < 4) that no other 
strong-disorder fixed point exist, the asymptotic behaviour is 
exactly given by the 1 -loop FRG equation (in d — 4) and by 
perturbation theory in the renormalized disorder, for d > A, 
i.e: 

A(u) = AoH-^/2[(Ao(u)-Ao(0))2]" + 0(A3) (B20) 



with h ~ {el2){I^ - 'm'^~'^)/{d - 4) (setting from now 
on K = 1 for simplicity). This equation is valid if we choose 
a model with weak and rough bare disorder, i.e. whose cor- 
relator exhibits a cusp. Of course, if we choose smooth bare 
disorder this perturbation formula fails, and one must run the 
RG to determine if the system is in the weak-smooth disorder 
phase where dimensional reduction holds (possibly up to rare 
events) or in the non-analytic phase. Finally, for d > 6 the 
situation is qualitatively similar up to additional dependence 
in the UV cutoff. 



APPENDIX C: GRAPHICAL INTERPRETATION OF K.V[w] 
AT 1-LOOP ORDER 



In this appendix, we give an intuitive derivation of the self- 



consistent equation ( 151 1 at 1-loop order. Two classes of dia- 
grams contribute, Ci and C2. We use the diagrammatics ex- 
plained in Section |V C| which can be used both for statics 
and dynamics. Here we chose to show explicity the arrows 
of causality. At this order, no difference is expected between 
static and dynamic shocks. We do not distinguish here be- 
tween Z and its rescaled version Z and use loose notations. 



These are the diagrams in the calculation of Z which, be- 
fore expansion in w, are proportional to A(w) — A(0). There- 
fore in order to get the term of order w, A(k;) — A(0) has to 
be expanded, and all other disorder vertices will be taken in 
the limit of w — > 0. For the n-th cumulant, we obtain (with 
A(w) - A(0) sitting at the bottom, and A" at the top): 



/C 



= nl 



1 



n-2=/i+i2 



-2A'H 



-2A'H 



/l>0,/2>0 

X [A"(u;) + A"(0)] [AH-A(O)] 
= n\ (n - 1) 2"-3[-A'(w)]"~2 [A" {yj) + A"(0)] 

X [A(u;) - A(0)] 
= -n\ {n - 1) 2"-2|A'(0+)|"-iA"(0)u; . (CI) 

li and I2 denote the number of A' in the left and right leg 
respectively. Rescaling by Sm amounts to set |A'(0^)| — > 1. 
Hence we find that the coefficient An defined in the text reads, 
for this class: 



(C2) 



Since in the self-consistent equation each outgoing line in the 
above is branched to a Z, this implies that the contribution 
of this class of diagrams to the r.h.s. of the self-consistent 
equation can be summed into: 



E 

n=2 



Ai^ 



I,.Z'' 



1 



-2ZY 



(C3) 



In the end one sets w? — ^ 1 for rescaled quantities. Starting 
the sum instead at n = 3, as we did in the main text, yields 
the additional subtraction —Z"^ J^, pjjqi^^,^)^ ■ This procedure 
identifies the latter term as counter-term of the disorder renor- 
malization, thus is the correct expansion in terms of the renor- 
malized disorder, whereas ( |C3| l is the correct result in terms 
of the bare disorder. 

We can now give a graphical interpretation to the dressed 



propagator l/(fc^ 



2Z) which appears here. In the 



above formulae, the left and right legs are effectively dressed 
propagators, which we represent as a double line: 



(C4) 



Class d: 



28 



The notation on the r.h.s. of the equation is as follows: The 
left vertex of each disorder is at 0, the right one at w (the 
choice being dictated so as to have a non-vanishing vertex un- 
der application of the K, operator). Aitows pointing to the 
right or left are going out - they will be branched to an exter- 
nal point sitting at (left) or w (right). From each disorder 
vertex A'{w), there are 2 outgoing lines, resulting in a com- 
binatorial factor of 2 for continuing either with the left (0) or 
right (w) vertex. Also note that the object exists, indepen- 
dently of whether the in- and out-going lines are at or w. 
Now class Ci can be written as; 



Ci = 



(C5) 



(Actually, the graphical notation is a little sloppy, since the 
upper vertex can either be A"{w) or A"(0), but there is an 




additional combinatorial factor of i in eq. (Cli, first line) 



The result is in agreement with ( C3 1 



Class C2: 

This class of diagrams looks like a correction to the critical 
force 



Co = 




= [A'H-A'(0+)] / 

Jk 



Z 



- 2Z) B 



(C6) 



Indeed it should be viewed as a circle of A', of which exactly 
one is expanded in w, leading to a circle with one marked 
vertex. This is the vertex drawn above. Note that the double 
line needs at least one A'(w) otherwise it cannot start at 
and go to w as indicated. This leads to the last term in (C6 1 
being subtracted. If one wants the expression in terms of the 
renormalized disorder, one again has to subtract the contri- 
bution proportional to ^^2^"^„^2-)2 , giving an additional term 

-[A'H-A'(0+)]/,(,4|^. 

(|C5[ l and \C6) (with the proper subtraction) together give 
( fl49r 



APPENDIX D: MOMENTS FROM Z{\) 



A direct series expansion of the formula (166 1 in powers 
of A, e.g. using mathematica, yields the following formula for 



the moments {S expressed in units of Sm) 

{S") (-4)"-iy7F 



{S) 



3(-l)"4"-2^i/' (i -n){n~ 1) 
(-l)»22»"37^(n- 1) 
3(-l)"4"-^V^(-2 + 7 + log(4))(n - 1) 

r(i-n) 

(-l)"22"-'''y^^ (i - n) 



^r(n) - 



(-1)"22 



r (i - n) 
(-l)»22"-30F(-2 + 7 + log(4)) 



which yields the universal ratios: 

/ (-i)"r(f-«)r(n-i) 
^ a\^n-\ — — 

7'n = 1 



(Dl) 



2n-3 



(2n-3)2 



(D2) 



We note that these formula make sense only for n integer, 
and in fact one can check that for integer n ^ 1,2,.. they give 
the same result as the formula in the text obtained from P{S). 
Because of the factor (—1)" they are not suited for analytical 
continuation to any real n. However using reflection identities 
of r(a;) and = ^'{^)/^{^) functions one can eliminate 
the factors (—1)" and get back the formula of the text which 
are real for all real n. 



APPENDIX E: CALCULATIONS FOR NON-LOCAL 
ELASTICITY 



Let us first recover the result for local elasticity 7 = 2, 
dc — 4: with 3^7^ = fc^ + 1. One can rewrite the momentum 
integral defined in the text as: 

Z2 Z 



J2 = N2 r dss"'^-^ 

Jo 



The constant N2 is defined as 

N2 = ^ 



(.s+l-2Z)2 .S+1-2Z 



z z^ 

5+1 (S+1)2 



(El) 



^ Jo "-^^ 



Explicit calculation yields back ( 149 1 and for d ^ 4: 



J2 = ^Z{2Z + (1 - 3Z) log(l - 2Z)) + 0(e) 



(E2) 
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i.e. one recovers ( 150 1 by a different method, which is also a 
check of ( |234l ). 

Let us now study the case of general 7 using the model 
(12281). 



1 



(1^1-^ + 1-2^)2 |fc|7 + i-2Z 



Introducing iVf 
computed as: 

Jo 



|fc|^ + l (|fc|^ + l)2 

2^(rd..''/^-i(^)-\ this 
^2 ^ 



(E3) 



(.s+l-2Z)2 S+1-2Z 



S+1 Gs+1)2 



1 ~ 



7 



Z(2Z + (1 - 3Z) log(l - 2Z)) + 0(27 - d) 



(E4) 



i.e. a function identical to (E2i at the critical dimension up 



to the global multiplication by 2/7. Hence the distribution 
of avalanche sizes to one loop with this choice of non-local 



elasticity will be exactly given by (681 with the replacement 
of a 2q:/7 in ( 170 1 and ( 171 1 as detailed in the text. 

Let us now study the case 7=1 and the form ( |229| l suitable 
to describe the contact line. We need to compute: 



{VWTI-2ZY 
z z 



(E5) 



Z2 



where = (e 



VFTT - 2Z ^/W+l k'^ + 1 

00 /c''"^dfcN~l 



fe2 + l 

following integral representation 



, € = 2 — d. We will use the 



1 



Vl + fc2 - 2Z 
This yields: 



00 -(l+fc^)Qt^+2Zi- 



0=0 Jt=0 



20Fa3/2 



(E6) 



1 



1 



VI + fc2 - 2Z VT 
r(i-f)r(|-i)bFi(i, 



2Z 1 - AZ' 



i-rf . 1. 
2 ' 2 ' 



2-d 
4^2) - l]sin(^f ) 



27r3/2 

Taking a derivative w.r.t. Z yields 



(E7) 



l-AZ^ {4(d-l)Z^-l 
{Vm?-2Z)^~' d-2 
4Zr(| - |)r(f - 1) 2i^i(2, 1 - I; |;4Z2)sin (f ) 



r3/2 



(E8) 



Together this gives 



Z^ 4(d+ 1)Z2 -3 1 -4Z 



3Z2 



zr(i-f) 



d-2 d-2 
2Fi(2,i-f;i;4Z2)-l" 



2Z^ 
1-2Z 



2VnT{2-l) 
- 3^2 log(l - 2Z) + 0{2 - d) 



(E9) 



(ElO) 



This result (ElO 1 has to be compared to (E4 1 for 7=1. One 
sees that the leading behavior for Z — > —00, from which is 
extracted the exponent r, is identical hence the avalanche size 



exponent r is still given by (237 1. The form of the avalanche 



distribution however will be different. To see this, first invert 



Z = \ 



2Z^ 



1-2Z 



-3Z^log(l -2Z) 



(Ell) 



This gives 

z!^^,W^\[i-Vi^) 

{VI^ - 1)2 [f Vr^log(l-4A) + - 1] 



4(1 -4A) 



(E12) 



One can now verify that to first order in a one can write (set- 
ting Sm ^ 1, restored in the text): 

/•oo 

Z^^i(A)=/ dS [e^'' -1]p!;^,{S) (E13) 
Jo 



with 

X exp 



1 



2V^ 
1 



1 + -(8-37e) 



5- 



^ l-^(7E-2)a)5i-^ 



(E14) 



which is the form given in the text. Note that to first order in 
a it can equally well be written as 



p!;=2{s) 



1 



2V^ 
X exp 



1 + -(8-37e) 



„_3_ 3a 
5 2 4 



Ul-lu-2)a]S'-- 



'ira „3 



(E15) 



This would suggest an abrupt jump in the effective stretched 
exponential decay exponent 5 to the value 3/2. However since 
the above results are valid for fixed S at first order in a one 
cannot conclude on this point at this stage. 
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APPENDIX F: MORE ON SPATIAL STRUCTURE OF 
AVALANCHE DISTRIBUTIONS 

1. Derivation of the recursion relation 

We want to show in Fourier space, using (t){x) — cos{qx): 



Zk{\) = A- [Sk+q + Sk-q] 



+ |A'(0+)| / gpgk-pZp{X)Zk-p{\) , (Fl) 
Jp 

where Sg = (27r)'*5''(q). 

The calculation of G"^(A) is a generalization of the one of 
Section [V] to non-zero momentum graphs. One still takes w 
uniform and now computes the spatially dependent correla- 
tions (for Qi all non-vanishing): 



9q^Ug{wi)Ug,{w2) = ~Rq[wi - W2]Sq+q' (F2) 



Uq, (Wl) ■■■Uq^ [Wn) = (-1)"^,'^ [wi] ' ' ' V^'^ [w„ 
= {-irCl^lqSw^.....W^)5q,+...q^ , 



(F3) 



where 5q = {2tt) 5 {q). We use Fourier space versions of the 
definitions in Q. In particular we have defined 



J X 



SWxSWy 



using translational invariance for a uniform w and the defini- 

— c 

tion of the R[w] functional 1621 y[u'i]y[w2] = R[wi — W2]. 
The q = limit is recovered using Uq=o{w) — L'^u{w). From 
these C one computes the associated momentum-dependent 
Kolmogorov moments, applying the operator JC. The choice 
(j>{x) = cos{qx) amounts to choose qi = ±q on the external 
legs of the tree diagrams of Section |Vj as for the calculation 
of IS'qP", hence we need qi — ±q. 

We extend the previous tree-level analysis and associate 
0"m,m' to the box with m times +q and m' times —q enter- 
ing from the top, and a moment (m ~ ■m')q exiting from the 
bottom of the tree diagram. Pasting together the trees as ex- 
plained in Section|V]we obtain the recursion relation; 



Em! m ! 
'^p^p''^t.i' p^\pnin9ip-p')q9{i-i')q 

p-\-l—7n; p' —m' 

(F5) 

Introducing bm,m' ■= cim.m' /mlm'l one obtains 
6™,™'e^^("-™')« = / g^gy (F6) 

Jxy 

p-\-l—m\p'-\-l'—m' 

By definition: 



In the sum ( |F5| l the terms 1 , and 0, 1 are left out, so one finds 
in real space: 

Z^(A) = Acos(gz) + |A'(0+)| / 5^ {X)Zy{X) 

J xy 

(F8) 



The Fourier transform of this relation gives (Fl 1 and it is a 
particular case of ( |202[ ). The general case is obtained by at- 
taching (t>{x) to each external leg. At the end of the calculation 
one wants Zk=o{^) since no momentum flows from the lower 



vertex, hence formula (201 1. Note that we have neglected ev- 
erywhere the non-local component (i.e. the momentum depen- 
dence) of the A'(0) vertex which we have taken at zero mo- 
mentum. Near c? = 4 this should be sufficient, as discussed 
again below. 



2. Generalized cusp-moment relation and check on the second 

order 



Let us start from the general relation: 

{u^{wi) - Wi){Uy{w2) - W2f = - / gxz9yz'R!xy[wi - W2] 

J zz' 

(F9) 



Consider now the expansion to second order A^ of (198 1. It 
yields: 



(F4) G^A)|; 



1 



Xe '^(■^) J xyzz' 



4>{x)m9xzgyAKAwyKz'm 



Taking a derivative w.r.t. w we obtain: 

Z'^(^X)\x2 [ 



(FIO) 

x)(t){y)gxz9yz'R'zz't[^ = 0+] 

(Fll) 

Comparing with its definition, and using the relation (S"^) = 
(S) L^'^ J dx(p{x) yields the result for the second moment of 
the avalanche amplitude at different points in internal space: 

= gxz9yz'R7z't[^ = 0+] (F12) 



an exact relation which is the non-local generalization of ( 36 1 
To leading order in e the third derivative vertex is local i.e 
R'^'^nlw] = Szz'tR'" {w), hence we find simply: 



= 2L- 



or m (/-space 



{S) 

W) 



9xz9yzR"'iO+) 



(F13) 



'^^'^'^ =2\A'{0+)\L-'glSq+q, (F14) 



This is what is predicted by the self-consistent equation 

{{Sfr)/{S) = 2|A'(0+)|L-'^ / ( / gxyHy))' (F15) 



X Jy 



(F7) 



a simple generalization of ( 36 1, valid to lowest order in e (or 
in the locality expansion). 
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APPENDIX G: DISTRIBUTION OF THE 
CENTER-OF-MASS FLUCTUATIONS 

One defines the generating function for the probability dis- 
tribution of the center-of-mass fluctuations as 



(Gl) 



/(A) = L-'^lncxp(AL''u(0)) = L-'*lncxp(A / u^(0)) . 

The result \\35) implies that, to one-loop accuracy, 
A2A(0) f^l^ A'(0+)2PA2p 



/(A) 



2 m-* 

A' A(0) 
2 m4 
1 

+ 2 



In 1 



(G2) 

A^(0+)^A^ \ A^(0+)^A^ 



Note that the term p = 1 is contained in the renormalized 
A(0). Since we work to one-loop order, we can now compute 
the integral in the second term in d = 4. The calculation 
yields 



/(A) = + W5(Am-4|A'(0+ 



)l) 



1 



-,{b) = -((l + 6)2ln(l + 6) + (l-6)2ln(l-6)-3fe2) 
64 56 ^8 510 
~ ~48 ~ 240 ~ 672 ~ 1440 



0{b'') 



(G3) 



with 54 = e/2 = l/(87r2) in d = 4. We note that the scale 
Sm = m~'4|A'(0"*")| = 0{em^'^^'^) naturally appears. In- 
deed we can rewrite the result for the generating function near 
d ~ 4, using the rescaled correlator from appendix [B] as 

r \2 



exp 



' (mL) 



d r \2 



87r2 



(1 - 2Ci) 



A^ A(0) 
-Ta'(0+)2 

52) + .9(A) 



.9(A) 
• (G4) 



A is now a dimensionless number The total factor in the expo- 
nential is obtained correctly to both orders 0(l/e) and 0(1). 
In the last equation we have used the FRG fixed-point equa- 
tion for A(0). The term 32 is zero if one uses the FRG equa- 
tion to one loop. However to get the A^ term to 0(1) requires 

injecting the ratio ^^^^''^ to two loop. These ratio are given 
in the Appendix and it is easy from there to find the number 
92 for each universality class. 



The result (G4 1 seems to indicate 119511 that the tail of P{u) 
m|*), with 5 = 2 - e(l - 2Ci)/2 < 2, con- 
result (5 = 3. It remains to be understood 

0, 



decays as exp 
trary to the d 

whether 5 first decreases and then increases again for d 
or whether the functional form is different, and the assump- 
tion that the In can be re-exponentiated is unfounded. A larger 
tail (6 < 2) seems to be counter-intuitive, since all connected 
moments in ( |G3| l above the second are negative. 

The main result ( |G4| i can be interpreted as a nice example 
of the central-limit theorem. There are effectively mL'^ inde- 
pendent regions |96|. In each regi on the generating function 
is given by the same formula ( G4 1 setting the mL'^ factor to 
unity. The final distribution of u(0) is obtained as a convo- 
lution of these mL'^ identical ones. In each region one can 
write u{0) — ^J~lu\ + eu2 — Smi^ui + U2) where mi,U2 
(resp. Ml, U2) are independent random variables of order one. 
Ml (resp. ui) has a gaussian distribution, and U2 (resp. 112) has 
a non trivial distribution encoded in the function .9(A). Hence 
near d = 4 the center-of-mass fluctuations are gaussian with 
a small correction encoded in g{X). 
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Note however that even if the weight in probability of 
avalanches of the size of the UV cutoff is negligible as m ^ 0, 
these may still control some (negative) moments of the proba- 
bility distribution. 

It depends on it subdominantly, i.e. in (5'min/5'm)'^^^. 
We use that j^^Au/\"{u) = 0, to 

conclude that — limj^o /.^^ dit A"(ii) = 

lim.^o [HL ^"(") + /r d« A"(u)j = -2A'(0+) 
We thank A. Fedorenko for this observation 
[93] To be precise, this means in powers of the local part R{u) of the 
functional R[u\, since higher multilocal parts can themselves be 
expressed as a function of the local part R{u), see e.g. 1631 . 
[94] The case of a vanishing J^(f){x) has to be considered separately. 
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from —1 to — oo, and from 1 to oo. Taking u large, only the 
left one matters. At leading order in e, one can still follow al- 
most the same contour, except that the integral from —u/2 to 
— 1 below the cut and from —1 to —it/ 2 above the cut has to be 
added. Its non-vanishing part comes from the 2ni discontinuity 
of the In across the cut, thus does not give a In contribution. The 
final result at leading order is then obtained by replacing A by 
A — u/2, doing the (almost) Gaussian integral in A in imaginary 
direction, with the result that P(it) ~ (G4i|a^-u/2- Exponen- 
tiating the In leads to the quoted result for 5. 
[96] Note that the above result is for mL ^ 1, the result for arbitrary 
mL can in principle be obtained by modifying the integrals into 
periodic sums. The opposite limit mL was studied in 1531 . 
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